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FOREWORD 


This  paper  describes  an  algorithm  developed  for  solving  a  special  class 
of  generalized  transportation* type  problems  of  moderate  size.  Problems  con¬ 
cerned  with  optimal  allocations  of  resources  subject  to  meeting  a  given  set  of 
requirements  such  as  marketing,  routing,  production,  and  weapons  allocation 
are  frequently  of  the  generalized  transportation  type. 

The  generalized  transportation- type  problem  considered  here  is  a  linear 
programming  problem  with  solutions  giving  the  allocations  of  the  jth  re¬ 
source  to  the  ith  operation  such  as  to  maximize  a  given  profit  function.  The 
requirements  specify  the  limits  on  each  of  the  n  available  resources  as  well 
as  the  operational  limits  of  each  of  the  m  operations.  In  addition  the  opera¬ 
tional  capacity  of  the  ith  operation  when  the  jth  resource  is  assigned  to  it  is 
known.  The  structure  of  such  problems  (one  constraint  for  each  row  i  and 
each  column  j )  enables  an  algorithm  more  efficient  than  the  general  simplex 
algorithm  to  be  used  for  finding  a  solution. 

The  algorithm  is  intended  to  solve  moderate- sized  problems  faster  than 
will  general  simplex  algorithms.  It  requires  less  computer  storage  than  gen¬ 
eral  simplex  algorithms,  thus  making  it  particularly  useful  when  a  limited- 
capacity  computer  memory  is  all  that  is  available. 


Nicholot  M.  Smith 
Head,  Advanced  Research  Department 
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ABSTRACT 


The  algorithm  described  in  this  paper  is  used  to  solve  a  special  class  of  linear 
programming  problems  characterized  by  constraint  coefficient  matrices  having  gen¬ 
eralized  transportation  structure.  Specifically,  n  available  resources  are  allocated  to 
m  capacity-limited  operations  (where  the  operational  capability  of  assigning  the  jth  re¬ 
source  to  the  ith  operation  is  known)  such  as  to  maximize  the  total  profit  for  the  system. 
The  row-and-column  structure  of  such  problems  permits  an  algorithm  more  efficient 
than  the  general  simplex  algorithm  to  be  used  to  solve  moderate -sized  problems  (prob¬ 
lems  where  loop-tracing  techniques  or  equivalent  schemes  are  not  required).  It  is  not 
required  In  the  problem  statement  that  all  the  resources  be  allocated  or  that  all  opera¬ 
tions  be  performed  to  capacity  limits.  It  is  characteristic  of  such  problems,  however, 
that  the  optimizing  solution  usually  requires  that  at  least  one  of  the  two  conditions  holds, 
i.e.,  either  supply  or  demand  is  exhausted.  The  paper  contains  a  description  of  the 
algorithm,  a  computer  program,  an  example  illustrating  its  application,  and  some  com¬ 
parisons  with  the  general  simplex  algorithm  in  solving  the  same  problem. 
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1.  INTRODUCTION 


The  algorithm  presented  here  yields  optimal  solutions  to  a  special  class 
of  linear  programming  problems  that  are  characterized  by  constraint 
coefficient  matrices  having  generalized  transportation  structure.  t  The  algo¬ 
rithm  preserves  primal  feasibility  and  the  complementary  slackness  condition 
at  all  times;  hence  feasibility  of  the  dual  constraints  forms  a  set  of  necessary 
and  sufficient  conditions  for  testing  optimality. 

The  need  for  the  present  algorithm  arose  initially  in  application  to  an 
optimal  weapons- allocation  problem  as  part  of  a  larger  nonlinear  minimax 
problem  employed  in  an  earlier  RAC  study1  in  this  area. 

A  specialized  algorithm  (similar  to  the  one  given  here)  for  generalized 
transportation- type  problems  appears  to  have  been  first  used  by  Ferguson  and 
Dantzig.2’3 

The  algorithm  can  be  divided  computationally  into  two  phases:  (1)  the 
matrix  maximum  phase  and  (2)  the  simplex  phase.  In  phase  1  the  algorithm 
permits  only  vectors  associated  with  positive  cost  to  enter  the  basis  and  only 
basis  vectors  associated  with  slack  variables  to  leave  the  basis.  In  phase  2 
the  selection  of  the  next  neighboring  vertex  is  currently  made  as  it  is  done  in 
most  simplex  algorithms  (see  Ref  4,  Lecture  V  and  the  appendix). 

The  particular  structure  of  the  constraint  coefficient  matrix  permits 
economy  of  computation  by  employing  the  equivalent  of  a  doubly  indexed  sim¬ 
plex  algorithm. 

2.  PROBLEM  STATEMENT 

The  algorithm  presented  in  subsequent  sections  yields  an  optimal  solution 
to  the  following  class  of  linear  programming  problems. 

tThe details  of  this  structure  will  be  considered  in  Sec  3,  ‘Problem  Structure.” 
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Maximize 


1  c'x'  with  respect  to  x‘ 


subject  to  the  constraints 


m) 


1  hrx'f  5  bj  ,  bj  >  0  ;  (j  =  I . n) 


i=l 


x  >  o 

»J  -  u 

2 0 

v.i>0 


(i  1 . m) 

()  =  1 . n) 


Under  the  correspondences 


xn  =  V-'. 

d .  =  r  /<r 

•I  ii  n 

cn  “  fu% 


(1) 


(2) 


an  optimal  solution  to  Prob  1  can  be  found  from  finding  an  optimal  solution  to 
Prob  3. 

Maximize 


m  ,rt 

i  c  ,x  with  respect  to 

i,/  "  •'  ’ 

subject  to  the  constraints 

■?,  V«i  -  ’  *i  >  0  ;  (l  =  1 . m) 

(3) 

m 

,S,  *„  <  bj  ,  b,  >  o  ;  (j  1 . n) 


The  algorithm  finds  an  optimal  solution  to  Prob  3. 

It  should  be  observed  that  both  the  row  and  column  constraints  are  inequali¬ 
ties.  It  is  characteristic  of  such  problems  that  the  optimizing  solution  has  the 
property  that  either  all  the  row  constraints,  or  all  the  column  constraints,  or 
both  all  row  and  all  column  constraints  are  binding  when  all  c()  >0.  If  equalities 
are  imposed  on  the  column  constraints  and  the  row  inequalities  are  of  either 
type,  we  nave  the  generalized  transportation  problem  considered  by  Hadley.5 
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If  both  row  and  column  constraints  are  equalities,  E  a,  =  E  by,  and  d,y  =  1  for 
all  i,j,  the  problem  reduces  to  the  standard  transportation  problem. 


3.  PROBLEM  STRUCTURE:  GENERAL  DISCUSSION 


The  general  simplex  algorithm  may  be  used  to  solve  Probs  1  or  3.  For 
large  m  and  n,  however,  it  is  not  practical  to  do  so.  Writing  the  components 
x(.  of  t  e  Emn  using  a  single  component  subscript  index  fe  for  xfe  (as  is  done 
when  using  the  general  simplex  algorithm),  we  see  that  the  constraint  matrix 
A  contains  mn(m  +  n  -  2)  zeros. 

If  g  is  the  component  subscript  indexing  function  [gO,  i)  =  fe]  for  the 
vector  x,  then  for  the  problem 

max  (c.x)  x  e  £mn  (  e  Em" 

A,  <  b  be  Em,n  b  >  0  (4) 

x  >  0 


A  assumes  either  of  the  two  structures 


(5) 


when  0(i,j)  =  n(i  -  l)  ♦  j  =  k  (1  <  i  ±  m,  1  5  j  <  n) 
where  In  is  the  identity  matrix  of  order  n  or 


(6) 


when  0(i,j)  =  m(j  -  1)  +  i  fc(i<jin,l<i5m) 
where  lm  is  a  row  vector  of  m  ones. 

tThe  right-hand  side  column  vector  b  here  includes  m  +  n  components  (m  a,  and  n  b() 
as  in  Prob  3. 
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Any  other  indexing  by  g  (besides  interchanging  upper  and  lower  blocks) 
produces  a  less  uniform  structure  for  A.  Structure  5  for  A  is  associated  with 
generalized  transportation-type  problems.  When  all  =  1  the  structure  of  A 
in  Structure  5  is  that  of  coefficient  matrices  associated  with  transportation 
problems.  Structure  5  for  A  will  be  assumed  when  introducing  suitable  basis 
vectors  for  the  solution  space  later  on. 


4.  OPTIMALITY  CRITERIA 


The  algorithm  preserves  primal  feasibility  Pand  the  complementary 
slackness  condition  S  at  all  times  and  uses  the  feasibility  of  the  dual¬ 
programming  problem  constraints  D  as  the  optimality  test  criterion. 

The  three  sets  of  Conditions  P,  D,  and  S  are  explicitly 


Ax 

+ 

1  i  =  b 

/A  w  -  l_„u-  = 

m+n  s 

1  mn  s 

X 

> 

0 

*  0 

> 

0 

IV 

o 

S:  <  ums  +  <u;s,i)  *  0 

where  '  denotes  transposition 
<  ,  )  denotes  inner  product 
b  :  0 
c  >0 

See  Ref  6,  Pt  2,  p  58,  for  a  discussion  of  Conditions  7. 

Real  vectors  x,w  that  satisfy  Conditions  P,  D,  and  S  also  solve  the  pair  of 
dual  linear  programming  problems 

max  (c.l)  subject  to  P  (8a) 

min  (iv,b)  subjecttoD  (8b) 

tv 

Problem  8a  is  solved  with  Conditions  P  and  S  always  holding,  hence  Condition  D 
becomes  the  set  of  necessary  and  sufficient  conditions  for  optimality. 

In  practice  the  algorithm  enforces  the  following  stronger  form  S  of  Con¬ 
dition  S,  namely, 

(w,xs)  +  (ws>x')  =  0  at  the  component  level.  (7S) 
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i.e.,  ijj>Oa»s^O,is  >  o  »  iVj  >  0  .  Since  the  dual  space  constraint 
A w  -  lmnu>s  =  c  holds  for  all  throughout  the  algorithm,  Conditions  7D  be¬ 
come  the  set  of  necessary  and  sufficient  conditions  for  optimality,  i.e., 

Ws  >  o  .  >  o  (7D) 

5.  PROCESSES  OF  THE  ALGORITHM 

The  general  processes  of  the  algorithm,  the  details  of  which  will  follow,  are 

(1)  Generate  basic  primal  feasible  solution  using  complete  or  partial 
matrix  maximum  (Conditions  7  P  are  satisfied  by  exactly  m  +  n  positive  primal 
variables  x,xs  while  the  primal  objective  function  is  increased). 

(2)  Solve  for  the  dual  space  variables  u,  and  ^  using  Condition  7S  and 
knowledge  of  the  structure  of  basis  vectors  associated  with  the  positive  primal 
variables. 

(3)  Perform  optimality  test  (test  Conditions  7D).  If  there  are  no  viola¬ 
tions  the  current  basic  primal  solution  is  optimal. 

(4)  For  nonoptimal  solutions  find  the  largest  violation  of  Conditions  7D. 
Identify  the  associated  vector  for  entry  into  the  basis  for  the  primal  solution 
space. 

(5)  Find  the  representation  of  the  entering  vector  in  terms  of  vectors  in 
current  basis. 

(6)  Preserving  primal  feasibility  Conditions  I>  select  vector  to  leave  the 
current  basis. 

(7)  Express  the  solution  in  terms  of  the  new  basis. 

(8)  Return  to  step  2. 

6.  DETAILS  OF  THE  ALGORITHM 

The  detailed  description  of  the  steps  of  the  algorithm  is  presented  here. 

1.  Generation  of  Basic  Primal  Feasible  Solutions 

The  matrix  maximum  method  of  generating  solutions  x(J  is  a  process  that 
makes  allocations  (assigns  values  to  x;. )  to  payoff  elements  c(J  of  a  matrix  P 
of  payoffs  as  follows: 
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Let 


where 


/•  kif|«iri.  -  0  uni!  *„■  .  0)  (, 


. it) 

(/  1 . m) 

ti '  II  *  ) 

iti '  in  .  1 


i  d  i 

i  1  "  " 


v,  h, 


(9) 


x|f).  and  xm  -j  represent  the  residual  “slack”  in  row  i  and  column  j  (Prob 
3)  after  a  set  {xl(]  of  allocations  has  been  made.  In  Conditions  7  xs  = 

<*lit' . *mn'  ■  Vl . V/  •  If  either  *•«'  0r  V,  =  0  then  no  further 

allocations  involving  row  i  or  column  j  can  be  made  since  either  the  ith  capac¬ 
ity  has  been  achieved  or  the  jth  resource  exhausted.  Initially  P  is  the  matrix 
of  all  ct/  since  x|r).  =  a,  and  xm  (  =  ^ . 

For  each  allocation  x,,  of  the  matrix  maximum  method  let 

ur  7;*  e  Fl  (10) 


then  choose 


ur 


1 


I  I 


akf 


l  1 


(11) 


This  choice  for  the  value  assigned  to  x^p  eliminates  either  row  k  or  column  l 
from  the  matrix  P  of  payoffs  for  the  next  iteration.  The  new  values  x^n-, 
for  xfef|.  and  xm-f  are  found  as  follows: 


If 


Ur 


fed 

dkt 


then  (12) 

Un'  Uu'"  VUf  0 

Vp  Vp  -  UP 
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If 


then 


(13) 


‘fen'  ‘(ill'  <,JiPIm'f 


xm't  Vf  Xkf  0 


In  the  first  case  Pnew  P  lcfej  I  M-  °l  <j  I . n> .  In  the  second  case 


p  iclfi*;.f  oi  (. 


,  m) . 


Xl  • 

If  =  Vp  then  an  arbitrary  decision  is  made  to  perturb  x.  .  by  a  small 

dfcP 

amount  epsilon. 

The  matrix  maximum  procedure  can  be  terminated  in  either  of  two  ways, 
by  exhausting  the  matrix  P  of  payoffs  (complete  matrix  maximum)  or  by  assign¬ 
ing  a  fixed  number  (less  than  the  number  of  iterations  required  to  exhaust  the 
matrix  of  payoffs)  of  positive  allocations  xl(  to  be  made  (partial  matrix  maximum). 

Throughout  the  matrix  maximum  iterations  exactly  m  +  n  elements  of  the 
vector  X  of  allocations  (x|; ,  are  positive.  The  vector  X  satisfies  Con¬ 

ditions  7P  for  primal  feasibility. 

The  matrix  maximum  procedure  proceeds  from  a  vertex  of  the  solution 
space  to  a  neighboring  vertex  as  does  the  simplex  procedure,  but  specifically 
it  proceeds  to  the  vertex  that  has  one  less  positive  slack  component  and  one 
more  positive  nonslack  component  (i.e.,  component  having  positive  cost  ct])\ 
hence  the  former  is  more  efficient  using  a  per  iteration  comparison.  The 
matrix  maximum  procedure  is  not  sufficient,  however,  to  achieve  optimality 
in  general. 


2.  Solving  for  Dual  Space  Variables 

The  vector  X,  resulting  from  application  of  the  matrix  maximum  process 
(partial  or  complete),  is  a  candidate  optimizing  point  since  it  is  an  extreme 
point  (Ref  4,  p  58)  of  the  convex  set  K  of  points  (*,is ) '  satisfying  Conditions  7P 

Ax  .  In,„xs  b  ,  x  0.  is  0  (7p) 
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The  linearly  independent  set  (a  basis)  of  m  +  n  vectors  corresponding  to 
X°  =  (x°,r°) '  (the  subvector  of  positive  components  of  X)  is  defined  as  follows: 


If  >  0  then  <1,1?,  +  ■*  a  member  of  the  basis 


If  x,>  >  0  then 


If  x  ■,  >  0  then 


m  +  j 


is  a  member  of  the  basis 
is  a  member  of  the  basis 


(14) 


Recall  that  in  the  matrix  maximum  process  if  >  0  then  not  both  xin>  >  0 
and  >  0;  hence  if  d ^  is  a  basis  vector  then  not  both?*  and  ?m+j 

are  basis  vectors.  Conversely,  if  both  xin>  >  0  and  >  0  then  x^  =  0;  hence 
if  t  and  eV  ,  are  basis  vectors  then  d,,?,  +  ,  is  not.  The  set  of  m  +n  column 

I  m+j  ij  I  m+j 

vectors  selected  from  the  matrix  (A,Im+n  )  using  Definition  14  and  denoted  by  B 
(the  ordered  matrix  of  such  column  vectors)  is  thus  linearly  independent  and 
satisfies  the  condition 

B.V°  =  b 


Hence  X  is  an  extreme  point  of  K. 

Corresponding  to  X°  satisfying  the  equation  BX°  =  b  is  a  vector  w°  satis 
fying  the  equation  B'u;0  =  c°  where  c°  is  the  vector  of  costs  (payoffs) 
associated  with  X°.  If  the  m  +  n  components  of  ir  are  written  w  = 


<U1 . um-  v\ 

c'  is 


.  .  ,  vn)'  then  the  scalar  form  of  the  equation  B'tr  =  c  or  u/B  = 


II 

jro 

if 

*ij  >  0 

(m+n)  equation* 

=  0 

if 

**’>  0 

(15) 

Vj  =  0 

if 

V,  >  0 

Solutions  to  Eqs  15  satisfy  Condition  7S,  <ir,xs>  +  (irs,  x>  =  0.  Since  xin>  >  0 
implies  u,  =  0  and  xm-j  >  0  implies  ^  =  0  then  (u%  xs)  =  0.  Similarly,  if 
xj(  >  0  implies  ivs..  =  d,^  +  -  c,j  =  0  then  (it's,*)  =  0. 

3.  Optimality  Test  (Testing  Conditions  7D) 

The  set  of  necessary  and  sufficient  Conditions  7b  required  for  optimality 
of  X°  is  rewritten  here  for  reference. 


t?,  here  is  a  unit  column  vector  in  E",+n  . 

?,  =  (0,0 . 1,0, - 0)'. 

t — ith  component 
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(7D) 


A  'u  I  n 

mu  s 


CD) 


The  scalar  form  of  Conditions  7D  is 


u 

f  ,  1 

I  mn  equations 

•1 

ir . 

0  1 

mil  inequalities 

•I 

0  . 

0  1 

m  *  n  inequalities 

(16) 


Let  (i^* . u®  ,  .  .  .  ,  v®)'  be  the  solution  to  Eqs  15.  Since  Eq  16  must 

hold  for  optimality  we  must  have  uf  i  0,  u°>  0 ,  v°  0  Kum  .  m  «  n)  inequalities ] 
where 


d  U° 

Mi 


(17) 


If  ic®  >  o  ,  u®  >  o  ,  i>®  >  o  for  all  i,j  then  X®  is  optimal  and  the  algorithm  is 

•J  ] 

terminated.  If,  however,  ui®  <0  for  some  i,i  or  u®  <  0  or  r®  <  0  for  some 

s(l  1  I 

i  or  j ,  an  improvement  in  the  solution  X°  can  be  made. 


4.  Nonoptimal  Solutions;  Finding  the  Greatest  Violation  of 
the  Dual  Space  Constraints;  Identifying  the  Associated 
Vector  for  Entry  into  the  Primal  Solution  Space  Basis 


The  greatest  violation,  V  ,  of  the  dual  space  constraints  (Conditions  7D) 
is  simply 


V  -  min  [min  IM0 
•il 


01  min  lu>?  <  0|  ;  min  <  0| 


(18) 


Depending  on  which  of  the  above  three  bracketed  minimums  is  largest  in 
magnitude,  the  corresponding  vector  c nosen  to  enter  the  new  basis  is  one  of 
the  three  types  of  vectors  d|;^  *  ^  ,  or  . 

5.  Finding  the  Representation  of  the  New  Basis  Vector  in 
Terms  of  the  Current  Basis  Vectors 

Consider  the  three  cases  (a)  V  =  u,®  ,  (b)  V  =  u?  ,  (c)  V  =  v?  that  can 

•  if  1  ) 

result  from  Eq  18.  *  The  vector  equation  to  be  solved  for  a  singly  indexed  sys¬ 
tem  is 

•+  j-* 

\  "H  or  >k  B  \  (19) 

^The  bar  denotes  the  minimizing  index  or  indexes  in  Eq  18. 
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where  y*  is  the  vector  of  coordinates  of  Ak  relative  to  the  basis  of  column 


vectors  of  B. 


Corresponding  to  cases  a,  b,  or  c  the  following  vector  equation  is  solved 
for  yjj ,  the  m  +  n  components  of  the  entering  basis  vector. 


(a)  '  ?»t7 


(b)  t-  Is  y~‘^(d  T  ♦  ?  )  ♦  S  )•'*•? .  +  I  y'1..? 

i  [  ’>l  hi  "+i  i"  1  ,  "*) 


x  .  >  0 

•I 


>  0  V,  >  0 


Equation  20  leads  to  the  following  three  sets  of  scalar  equations  in  y]j,  yjj1’,  or 


mj  . 
y>) 

For  Eq  20a 


*  >!!<*„  +  if  1  =  1 
>0 


(i  =  I . m) 


1  y‘»  d 

j  '•»  0 

*U  >  0 


0  if  i  if  I 


S  y[{  ♦  yjj'j  =1  if  )  =  J 


x|(  >  0 


(j  =  1 . n) 


0  if  J  V  j 


*.(  > 0 


For  Eq  20b 


1  if  I  i 

V  vT"'  d  =  ■ 

(  7<f  i| 

x|(  >  0  lo  if  i  d  i 


(l  =  1 . m) 


For  Eq  20c 


\>0 


(j  1 . n) 


T>H,du  0 


(l  =  1 . m) 


*■1  "*  ° 


1  if  1  i 


0  if  I  t  j 


(j  1 . ri) 
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6.  Selecting  the  Vector  To  Leave  the  Current  Basis 


Once  the  vector  that  enters  the  new  basis  has  been  found,  the  associated 
positive  components  of  the  new  primal  solution  must  also  satisfy  Con¬ 
ditions  7P  for  primal  feasibility.  Hence  we  have 

(24) 
or 


(25) 


i)\e 

o 

\ 

0  ,0 
nr  w 

(a) 

d--c-  <  c 

i  |  i  m  *  | 

(n)  <f;J?7 

1  — 

B.\ 

0 

0 

(b) 

-* 

fT 

.  0 

lb)  f7 

(c) 

| 

(c) 

-* 

fmt7 

(a) 

»> '' 

<“) 

-♦ 

c  - 

BX° 

0 

(b) 

ly"' 

4  0 

(b)  ?7 

(c) 

Bv"” 

(c> 

(a) 

Hy’f 

(a) 

U 

> 

0 

(b) 

0y7"' 

.  o 

(b)  t- 

(c) 

dym'~ 

(c) 

(26) 


(27) 


Since  X°ew  >  0  we  have  in  particular  (a)  iP-  =  6  >  0  ,  or  (b)  x? .  =  0  >  0,  or 
(c)  x®--  =  0  >  0  corresponding  to  the  new  basis  vector  a,  b,  or  c.  The  remain¬ 
ing  m  +  n  -  1  column  vectors  of  Bncw  are  determined  by  eliminating  that 
column  vector  of  B  whose  new  associated  primal  solution  component  x° 

' '  n  r  vv 

vanishes. 

This  elimination  is  accomplished  as  follows.  Writing  the  expressions  in 
the  left  braces  of  Eq  27  in  component  form,  we  have 


X 


o 

•J 


new 


(a)  0y\\ 

<b)  0y\* 

(c)  0ym'l 
7  1 1 


for  all  (i,j)  3  0 


(28a) 


fi-  *?„■  ">!,!'  for  all  (l.ll  )  3  I  -  0 


m»nrw  m> 

,  or  x 


tly'  I. 
m  ) 


for  nil  3  t  *  0 

m  j 


in  ,  or  X 

lirw 


m  (nr 


(28b) 

(28c) 


Since  we  want  x° 
from  Eq  29 
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to  vanish  we  select  positive  0 


\ 

■  1 1 

If  the  minimizing  indexes  in  Eq  29  are  (i  ,j)  =  < p, q)  then  x®  =  0  and 
<fp((c*p  +  tm  q  leaves  the  basis  if  p  ^  m'  or  q  t  n’ ,  ?p  leaves  the  basis  if  a  =  n', 
and  leaves  the  basis  if  p  =  m’ . 

7.  Expressing  the  Solution  in  Terms  of  the  New  Basis 

The  new  solution  X°cw  has  components  expressed  by  Eqs  28a  to  28c  with 
6  =  8  .  In  particular,  as  mentioned  before,  x®  =  0  and  x?-  =  8  for  the 

P‘<non  I) 

primal  variables  associated  with  the  leaving  and  entering  vectors  respectively. 

8.  Return  to  Step  2 

Self-explanatory. 


.  EFFECT  OF  NEW  SOLUTIONS  ON  VALUE  OF  OBJECTIVE  FUNCTION 


There  is  associated  with  any  violation  of  Conditions  7D  a  new  solution 
( X*J  to  Conditions  7P  that  improves  the  value  of  the  objective  function  (c°,  X°> 
and  at  the  same  time  eliminates  the  specific  violation  of  7D. 

Recall  from  step  6  of  the  algorithm  Conditions  7P  are  preserved  when  a 
new  vector  Afe  enters  the  basis,  thus 


II X® 

°\  ' 

b 

l)\® 

"»>fe  • 

b 

\  »>k 

(30) 

l)(\® 

b 

For  the  corresponding  expression  to  the  objective  function  value  we  have 

<r°.<X°  -  //»„>>  •  »c. 

or  o  o  1  (31) 

(c°,\°)-  (Ac  ^  )  (non  objective  fund  ion  value) 

The  term  (cfe  -  <c°, yfe ))  corresponds  to  (ch  -  zk)  in  general  simplex  notation 
and  in  the  notation  of  this  paper  to  (a)  -u’s  for  1  ?  fc  "  m»i  when  Afe  =  d^f,  .  ?mtJ 
fe  =  n (i  — 1)  +  j,  or  (b)  -u,  when  Afe  =  c(  mn  ■  k  .  mtor(c)  r  when  A);  eMtJ 
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mn+m<fesmn+m  +  n.  Thus  for  positive  0  and  any  violation  of  Conditions  7D, 
i.e.,  u?s  <  0  ,  u,-  <  0  ,  or  <  0,  there  is  an  associated  improvement  in  the 
objective  function  value  of  magnitude  (a)  -Bit's  ,  (b)  -0u,  ,  or  (c)  -80.  when  the 
vector  (a)  Ak  ,  (b)  ,  or  (c)  em+j  enters  the  new  basis.  Condition  7S  guarantees 

that  the  violation  will  be  eliminated  for  the  next  iteration. 

Throughout  the  algorithm  values  for  zfe  (zk  =  <c°,  yfe^)  are  not  computed 
using  the  yk  representation  of  ^k  (i.e.,  the  representation  relative  to  basis 
vectors  B),  but  from  the  relation  zk  =  (w®,  Ak)  which  makes  for  greater  effi¬ 
ciency  in  computation. 

8.  COMPUTATIONAL  EXPERIENCE 

The  algorithm  briefly  called  MATMAX  was  originally  used  to  solve  the 
linear  subproblems  described  in  Ref  1  with  m  =  3,  n  =  4.  During  the  process 
of  convergence  to  a  single  larger  nonlinear  programming  problem  solution  to 
which  the  linear  programming  Prob  3  is  only  a  constraint,  it  became  necessary 
to  solve  the  linear  problems  in  the  order  of  ten  thousand  times.  The  need  for 
an  algorithm  faster  than  the  standard  simplex  algorithm  thus  arose. 


TABLE  1 

Solution  Times  for  MATMAX  and  Standard  Simplex  Algorithms0'3 


Number  of 
constraints 

MATMAX,  sec 

Simplex,  sec 

Simplex,  ^  sec 

m 

n 

5 

4 

0.K) 

0.60 

0.36 

10 

12 

5.56 

31.51 

15.84 

18 

24 

24.88 

297.10 

114.98 

^Solution  limes  ore  based  on  single  precision  operations  in  FOR  TH  AN  1\ 
using  the  IBM  7044  computer. 

^See  App  C. 


The  MATMAX  algorithm  has  been  compared  for  solution  time  with  the 
simplex  algorithm7  and  an  even  faster  simplex  algorithm  given  in  App  C.  The 
A  matrix  (with  identity)  requires  19,988  =  42  x  474  storage  locations  for  the 
m  =  18,  n  =  24  simplex  algorithm,  thus  limiting  the  size  of  “incore”  compari¬ 
son  of  the  algorithms.  Solution  times  for  MATMAX  and  standard  simplex 
algor ith  .  are  shown  in  Table  1. 
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The  success  of  the  algorithm  currently  depends  on  being  able  to  solve 
the  m  +  n  linear  equations  Eqs  15  in  u,  and  i>(  sequentially,  i.e.,  on  solving  for 
the  nonzero  u,  and  v.  in  terms  of  zero  valued  u,  and/or  i>j . 

If,  however,  it  is  not  possible  to  solve  the  system  of  Eqs  15  sequentially 
during  some  iteration  of  the  algorithm, an  attempt  is  made  to  bypass  the  diffi¬ 
culty.  The  algorithm  then  attempts  to  proceed  to  the  optimum  by  avoiding  the 
particular  vertex  for  which  Eqs  15  could  not  be  solved.  Currently  the  algorithm 
returns  to  phase  1  (the  matrix  maximum  phase), this  time  assigning  one  less 
x.j  having  positive  cost  c|(  and  one  more  positive  slack  variable  (partial 
matrix  maximum)  than  was  assigned  the  prior  solution  of  phase  1.  Beginning 
with  this  solution  in  phase  2  (simplex)  a  new  sequence  of  vertexes  is  generated 
for  which  the  problem  of  nonsequential  solvability  of  Eqs  15  is  frequently  avoided. 

The  above  process  has  worked  successfully  on  most  problems  of  moderate 
size  but  has  failed  on  one  with  m  =  30,  n  =  32.  In  some  situations  more  than 
one  return  to  phase  1  may  be  required  in  order  to  find  a  sequence  of  vertexes 
for  which  Eqs  15  may  be  solved. 

It  is,  of  course,  possible  to  solve  Eqs  15  when  sequential  methods  fail. 
However,  the  logic  of  loop-tracing  techniques  required  in  such  situations  is 
complex  and  is  not  currently  employed.  Alternative  methods  that  use  the  se¬ 
quential  solvability  of  Eqs  15  as  a  secondary  criterion  fox  selecting  the  next 
neighboring  vertex  are  under  investigation. 
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Appendix  A 

COMPUTER  PROGRAM  FOR  THE  ALGORITHM 

INTRODUCTION 

A  listing  of  the  computer  program  for  the  algorithm  follows.  The  algorithm 
has  been  used  successfully  on  problems  of  moderate  size  (see  Sec  8,  “Compu¬ 
tational  Experience”).  The  comment  cards  identify  appropriate  subsections  of 
the  program  as  described  in  Sec  6,  “Details  of  the  Algorithm.” 

The  success  of  the  algorithm  currently  depends  on  being  able  to  solve  the 
m  +  n  linear  equations  Eqs  15  in  u,  and  e,  sequentially,  i.e.,  on  solving  for  the 
nonzero  u,  and  i>(  in  terms  of  zero  valued  u,  and/or  i’( ,  as  discussed  in  “Compu¬ 
tational  Experience.” 

Several  working  arrays  are  used  for  calculations;  namely,  the  AE,  ABASIS, 
YY,  ABAR,  BBAR,  U,  V,  ISUM,  and  JSUM  arrays. 

The  value  of  the  objective  function  is  printed  out  every  kth  iteration  by 
setting  WRITE  =  k  on  the  first  input  card.  In  addition,  a  detailed  printout  will 
be  given  every  kth  iteration  by  setting  ITAB  1  on  the  same  card. 

The  computer  program  of  App  C  accepts  exactly  the  same  input  cards  as 
the  following  program,  with  the  exception  as  stated  in  App  C. 

The  subroutine  TODAY  called  twice  in  the  program  is  used  for  timing 
purpose  only.  The  general  user  should  not  call  this  subroutine. 


PROGRAM 


FORTRAN  SOURCE  L  1ST 

1  bw  SbUKCt  STATEMENT 

A  I  rt  E  1  b  H  A  1  M  A  X 

1  JML  NblUl  LF  (  33,  35)  »Ofc  (  36,361 ,  At  1  36,361  ,  Ad  AS  I  S  (  36 ,  3  b )  ,  Y  Y  {  3b,  36  ) 

2  JMtNjIUn  Ad  Ar  (  36)  ,di3Ak(36),ISUM(3b),JSUM(36) 

3  OlMtNSlo  N  A  A  (  3b)  ,BBI  3b) 

4  JMENSlun  XX  (  3  5  ,3  3  )  ,  U  (  3  5  )  *  V I  351 

5  1S('.!1  FORMAT  (  1  HI  ) 
c  isj.2  Format  (in  ) 

7  1  Si"  3  FORMAT  I4I1i.) 
i:  IS  s'- 4  FUKMAT(SrU.b) 

11  lSv  •  c.  F|iK«AT(l5A,bcl6.6) 

12  19.  :t  FokMATI 15X,4HMS=  ,12,4X,4HNS  =  ,  I  2 , 1  5X  ,  3CHFK  ED .  UF  UBJ.  FUNCT.  PKIN 

1  Tool  , 12,  CX, 2 SH DETAILED  PRINTOUT  IF  1TAB  =  1. 5X, 5HI TAb= ,  12// I 

13  IS  /  FOrvMAl  (  3  ;X  ,46nl  NPUT  CONSTANTS  FCR  OPTIMAL  ALLOCATION  PROBLEM//  ) 

Is  lv.  H  Pui\MAT(i5X,24HlNPUT  VALUES  FUR  CEI1,J)//I 

If  Is.  S  FurMA f< 15X,2sHlNPUT  VALUES  FCR  DEII,J)//) 

lb  1).  1  _•  Fo.<rtAl(13X,45HlNPUT  VALUES  FOR  A  A  ( 1  )  OF  THE  ROW  CONSTRAINTS//) 

17  |S  11  F  jR.«lA  TI  1  3X  »4  3r(I  isPUT  VALUES  FOR  BB(J)  OF  THE  COLUMN  CONSTRAINTS//) 

l  19.12  F.JRMAK  S^HOF  Sfc  Nt  R  AC  Y  OCCURS  FOR  ABAR(I)  BBAR(J)  TkY  NEw  DtLTAl) 
cl  IS  lj  FoRisATl  /HlMAX  =  ,13,7HJMAX  =  ,13) 

22  lSd4  FORMAT!  I 3HAd  AR ( I RA  X )  =  , E 14 . d , 1 3 HB BAR  I JMAX )  =  ,El4.d) 

22  Is. IS  ForMAT! //37H TuTAL  EXECUTION  TIME  FCR  ALGORITHM  =  ,  F 12 . 7.  IX, 4HSEC . ) 
c  4  1  S  ■'  1  b  F  J  R  M.  A  1  (  B  1  1 1  ) 

c  f  1  St.  1  7  FurNATC4'  X.23HOETAILS  UF  THE  SCLUTICN//) 

2b  IS..  IH  F  ijrmAT  I //44HNUMBBK  UF  ITERATIONS  AFTER  INITIAL  SULUT  ION  ,1s) 
cl  lS'.ls  F  U  'Mrt  T  I  //34HVAI  i  OF  PRIMAL  OBJECTIVE  FUNCT  I  ON,  E2o  .  b ) 

3  w  IS;  2  1  FokiAI (//32HVA  l  OF  DUAL  OBJECTIVE  FUNCT ION, E2C .6 ) 

31  IS. 21  FORMAT! 1;X,4HR  ,  ,  I  2 , SX  ,  7HC CLU MN  ,  I  2 , 5X , 1 1 H ALLUC AT  I  ON  ,L12.6,5X,23 

IHrl  TUKii  FROM  1  .uCATIuN  ,bl2.o) 

3 c  IS  32  FORMAT < //74HV  tES  OF  THE  DUAL  SPACE  VARIABLES  (LAGRANGE  MULTIPLI 
l  F  rt  S  ,  SHAGO*.  kICESI) 

33  IS  2s  FoKMMd  X.4HRUR  ,  I  2  »5X  ,6HL(  I  I  =  ,E12.6) 

3s  IS  24  Ft.KMATC  7x,7hC<JLOMN  ,  I  2 , 5X  ,6HV  ( J  >  =  ,F12.6) 

;  f  15i25  rOKMAT! //37HMAXI MUM  VIOLATION  CF  CUAL  CONST  RAINTS, E 13.6  ) 

3  fc  1  S' ’2o  F  ut\ M  A  T  I  / / 1  u  X  ,2 tiHS  TAR T  I  NG  NEW  SECUENCE  NUMBER.  14) 

37  1  S ..  2  /  F  uR.MAT  I  3X  »4h  R  J  it ,  I  j  ) 

4.  IS, 2b  FOR3ATI1  'X, 41 H TOTAL  1TLRATIONS  IN  LAST  VERT  ICE  SEUUENLE,I3) 

41  IS  ,2b  FuRMAKl  X,4/HNLMhER  OF  POS.  XX(I,J)  TO  BE  ASSIGNED  BY  M  A  T  M  A  X ,  I  3  ) 

42  IS.  3b  F  iM.i  A  T  (|5^,  10H1  TER  AT  I  UN  ,  14, 8X,28H  PRIOR  VALUE  OF  UBJ  FUNCTION  ,E12. 

ib.rtX,  Jlt-MAX.  VIOLATION  CF  DUAL  CCNSTR.  .E12.N/I 

43  IS.  31  FlkmAT(//32HMATkIX  MAXIMLM  ITERATION  NUMBER  ,14) 

4s  IS. 32  FOrMAT! 5X, 1 9H AL LOC A T 1  UN  SELECTED, 4X, 4 HKOW  ,  I  4 , 4 X , 5HCUL .  ,14,4X,SHX 
1  X(  I  ,  J) =  ,  F 1 2 . b  ) 

4  f  is;  33  H  ;  j  t<  A  I  (  SXl  3HENTFKING  VECTOR, 213, BX, 14  HLEAVING  VECTOR, 2  15,  dX,  20HXX! 

1  I l  ,  T t K ,JF UTERI  =  ,Fl2.6//) 

4c  1  S  ■  j  3  4  i- ur  *  4  I  (  5x,3bHVALCE  OF  OBJECTIVE  FUNCTION  =  ,E12.6) 

47  1  S'.  35  FuRMAT  (  /  //JfiX.jiiHDt  IA1  LEO  INTERMEDIATE  PRINTOUT/) 

SO  IS. 3b  FukmAT (  / 1 (  X,3‘HCORRFNT  SCLUTICN  ARRAY  XX!I,J)/) 

si  1S..37  F  URMAT  (  /  i  t>  X ,  1.  SHUNUSEl)  K  E  SOUKC  E  S  ,  5  X  ,  26HCOL  JMNS  1  THRU  NS  IN  URUER/) 
f2  IS  3d  FORMAT! /15X,1 7HONUSEO  C  APAC  1  T  I  ES ,  5X  ,2  3HRUWS  1  THRU  MS  IN  ORDER/) 

53  IS  3S  fukMA  T  (  / 1 3X ,  5  .11  PRIOR  VALUES  OF  THE  DUAL  VARIABLES  U(I)  1  4  URDER) 

5  4  1 S  '  4  FuRMATI  IV  X,5  H  PkIUR  VALUtS  OF  THE  DUAL  VARIABLES  V(J)  IN  OrDLK) 

5  5  RKl  TE  I  <  ,  IS'  L’1  ) 

ft  Wi\ITcIb,lVvt!7) 

f  7  KL  AD  I  3 » l  s(  3)MS,NS,IwKlTE,ITAB 

6  4  RKI Tt 1 b, 1 Su0  6) MS, NS, I WRI  TE  ,ITAB 
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FURTRAN  SOURCE  L  1ST  MA 1  MAX 


1  Sv 

SUJnoc  STATtMkNT 

6  S 

*r.  1  T  k  C  6  ,19'.'  2) 

t 

«  k  1  1  L  I  (> ,  1  4l  3  ) 

6  7 

UO  10..  5.  1  =  1,  MS 

7 

R LAD (3,1  on  4 ) (CL  I  1  , J)  , J  =  1 ,NS ) 

7  5 

rtkl  U  (0,14.  2  7)  I 

/  6  )  Si  V 

rtkllt  (o,  10'/  '  5  )  t  C  E  (  1  ,  J )  ,  J  =  1  ,  NS  ) 

1,4  „Kl  TE(o,19  ■  2  I 

1  *  «ikI  It  <  6,  1  9t  0  9) 

1  c  JU  L  f  v-oi  I  =  1  ,  M  S 

1:1  *t  Al.(  5, 1  *(  '.A)  (  Dfc  (  I  ,  J)  ,  J  =  1  ,NS) 

1  i 4  WKl  T  l  (  6  *  1 4<  2  7  )  I 

lit  ISf  W*l  Tfc(t>,lS\’  'illUK  I  ,J>  ,J=l  ,NS) 

12:  w  k  I  f u  ( /  ,  1  4 .  21 

12  4  *»' r< I  T i  C  t  » l  s (  l  j ) 

lei  RLAol  5,14''.  4)  (  AA(I)  ,1  =  1  ,MS) 

i  j  2  rfkl  Ult,,  l«Ji  till  AA(  I)  ,1  =  1  ,MS) 

1  3  1  wk  I  TCI  6,  ISC  2) 

14','  TL(o,14i  11  ) 

141  KL  Au  ( 1. ,  l  O'  J)  ,  J=1  ,NS) 

1  4o  Wkl  T  L  (  o  *  1  5 '  .'5)  <Bt)<  J)  ,J=1  ,NSI 

1  6 i  ri!<l TL ( o, 14w  1  I 

it.4  call  id j*.yc:,i  nut, iuati 

C - OL  <ltKA  I  ll..l  Lk  BASIC  FEASIBLE  SOLUTION  USINO  MATRIX  MA  X  1  MUM-  17  l  Jv-  1  99 

155  m  Total  =  ms  *  NS 

lit  MHlu  =  .VS  ♦  NS 

15?  OlLTAI  =  .ll-i 

U  y  JL  L  1 4  2  =  .  1 1  -  •/ 

10  1  :i  S  S  =  U  S  ♦  l 

162  MSS  =  vs  ♦  1 

16.1  MOOT  =  M  b  1  (, 

16  4  LSI  ‘ «  (1=1 

1  c5  T  L  1 P  7  =  . 

166  MMaT  = 

It?  i)  J  1/1-5  I  =1  ,MS 

1  7 1  S  1  '  5  A  t  A  i<  (  I  )  =  A  A  I  i  1 
172  JO  l'UH  J  =  1  ,  N  S 

1  7  j  i  Si  1  ilhAK  1  J  )  =  0*  (  J  t 

175  |)U  1/lr)  L=1  ,  MU  C  T 

176  A  MAX  = 

177  %  1/14  1  =  1  *  M S 

2  O  iriAvH.I  1) -.16-6)1  014  J,  1014;  ,  1912./ 

2.1  1512.  UO  1  /1j5  J  =  l  , 'M S 

2J2  IHHnV.lJ)  -  .Ik -6)14155,191  15*  10125 

c,3  1012  5  II  I  C  c  I  1  ,  J  )  -  AMAX) 191 15 ,19135, 191 ll 
<C4  1S1  1(  A  6.  A  X  =  Ct  (  I  ,  J) 

2'  5  I  MAX  =  I 

2  .  t  JkUx  =  J 

evl  10135  Cl.  /T  I  MUl 
211  l4l4t.  Cu.mTIhLc 

eli  1F( A MAX) 192  J , 1 02 . o , 1 ? 1 5u 

2  14  1 0 1 5v  AoAnTP  =  AuAkIIMAXI/  UE < I  MAX , JMAX ) 
t15  Ac(  1,1AX,JMAX)  =  Ctl  IMAX  ,  JMAX  ) 

<16  Brink.  TP  =  BBAKIJMAX) 

21/  lFUrtARTP  -  HBAKTP)  19160  *  191  70 , 191b1’1 

C - BRANCH  191  7t  IS  FOR  DEGENERACY - 


20 


SuUkLC  ST  All  He  NT 


F  UK  T RAN  SOURCE  L1SI  MATMAX 


I  S  4 


22..  Isl6t.  XX(  1  MAX,  JN'AX)=ABAKTF 

2  2  1  AAAk  (  I  lAX  )  =  _  . 

del  (Li  An  (  JRA  X  )  =  URARTP  -  AhAPTP 

22i  IFtlTAd  -  1)  1/H9, 19165, 19199 

c24  1  5  1  e  5  TI/MP2  =  Tt  MP2  ♦  C  t  (  I  MAX  ,  JMAX  )  *  X  X  (  I  M  AX  ,  JM  A  X  ) 

2  2  5  1 A  T  M  A  =  MM  AT  ♦  1 

2  2  6  wK I  Tc  (  t  ,  i  c<  3  1  )  MA  TMA 

22  7  WnlTC(7,i8'.i?)lMAX,JiM4X«AXIlMAX,JMAX) 

2l _  WKi  TL<6,15<  34)  TEMP2 

dll  1,1  III  1  ,1  O', 

dU  1517c  A I  ( A  r,  (  IFiAX)—  AGAk(IMAX)  ♦  CtLT 41 
21 3  lL>tGE9=ll>CGcN*l 

dl<i  iKUtbt'J  -  MTCTAL)  1915C  ,1915u, 19175 

2  35  1S1  75  rtKl  TH6,19i  1?) 

dit  wkl  TM6,i4t  l  31  IMAX.JMAX 

2 11  rfftl  TL  <  C, 1  1  A ) AUAK ( 1  MAX  1 ,BHAR<  JMAX ) 

7  4  2  1  Lib  Gt  14  =  ( 

<-41  JFLlAl  =.  13.  *  OELTA1 

2  4  2  GuTu  1917c 

243  ISlfcC  XXI  1  MAX,  JMAX  1  =  '3BAKTP 
(  i A  HuAk(JMAX)  =  J . 

24s  AtUnllMAX)  =  AdAHUMAX)  -  Ofc  (  1  MAX  ,  JMAX  (  ♦  BRAkTP 

?4t  IFIITau  -  1) 19199, 191d5, 19199 

<r  4  7  1  5 1  H  S  TEMP2  -  TIMP2  ♦  CF ( IMAX  ,  JRAX)  *  XX  (  1MAX  ,  JMAX  ) 

2  5.  MAMA  =  MmAT  ♦  1 

25  1  wK  I  Tc  (  r, » 191  31  )  hATMA 

2  52  *kl  TEtt>,l‘.K  1 2  )i  MAX  ,  JMAX  , XX  ( I  MAX  ,  JMAX  ) 

25  3  rfK I Tc  t  6, 1  Sc  a4> TFMP2 

2  54  1  5199  MMnT  =  HAT  ♦  1 

C - SOLVE  6uF  I  Ht  UUAL  SPACE  VARIABLES  Jill  AND  V(J)  192*0-19295 

2  5 1  152'..  ILUAl  =  . 

c57  JO  192L6  1=1 ,MS 

2c.  IFUJAnU)  -  .lfc-6) 192c 2, 19292, 1922 4 

*6  1  152.  2  IJI  1  )  =  1  .  fc  ♦  3  ‘j 
262  50  T  J  1921 fc 

2  6  3  1  525  4  UU)  =  f  . 

2  04  I  F  (  AA  (  1  I  -  AliARl  1  ))  192i.'f  ,l92-'6, 192l  5 

2 1  5  1  5 2 1 5  I  DUAL  =  IDLAL  ♦  1 

26F  152ct  CdNTlilOL 

27.  )u  15212  J= 1 *NS 

2  71  IFlbdAxlJ)  -  .  lt-6)  1921.6,19208,  1921C 

27c  15c 6  V(J)  =  l.l»35> 

2  ?  5  GO  r.J  19212 

2  74  1  52  1  r  V (  J  1  = 

2  7  5  IFIV.IJI  -  G  ti  Ak  (  J )  I  19212*  ,921  2*  1921  1 

276  15211  luUAl  =  IjOAL  «■  1 

2/7  152U  CunTImjl 

3  1  It  (IDEAL  -  1  1  19215,1921  5,19219 

3 .2  15215  MliUl  =  MM  A  T  -  1 

;.  3  IFIMJOI  -  c  I  1941  1 ,19216, 19216 

2  s  15216  JU  19718  1=1, MSS 
5  >u  1971  l  J  =  1  ,  NS S 

vc  A I  I  I  ,  J )  =  u  . 

1.7  1  5c  1  7  C un T  1  xOE 
ill  152lo  lUNT  I  NLL 
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FUkTKAN  SGUkCt  LIST  MAT  MAX 


I  SU 

Si  nJKuL  SIAIlNENT 

J  1  b 

1  Vtk  I  =  l  Vt  K  T  ♦  I 

ia 

WKl  Ttln,  1  S'  2fcl  l  VCkT 

.*  lb 

m/hl  TtlkfiSv.  23111 

-  It 

Wkl  T  b  ( t ,  1  s ^ H G U T 

:•  1  i 

1  Pit  1  i  1  = 

Go  Til  191  ft 

hi 

1  v219 

l  M  N  A  L  =  l 

ill 

'I  UJAL  =  '7 

Hi 

1  S??t 

1 1  (  IF  INAL  -  MS)  1 922  2  ,193Su»193uO 

32A 

i  C.{1 21 

Ml) UAL  =  Ml.  UAL  +  1 

S 

KIMJUAL  -  MbIG)19225,19225,1921b 

1  i  t. 

1  S22b 

1  F  1  ix  A  L  =  ( 

■11 

■JlI  1  s  2  S '  1  =  1,  MS 

“  ^  . 

IK  U(  l  )  -  1  .F+35) 19260,19230, 1926* 

bbl 

1  S  2  3  ■' 

l)U  1  9i_  S  J=1,MS 

332 

IK  V(J)  -  1  .E  +  35)  1923b  ,19250, 19235 

3  3  3 

1  S?  jb 

If ( Atl I , J) ) 19Ebs  ,  1925o  ,19240 

32A 

1  S2  A  l 

'J(  1  )  =  (C.t  (  1  ,  J)  -  V  (  J)  )  /  DEI  I ,  J  > 

O  i  c; 

GU  Tu  1  9  2  S  i 

'.J  It 

192  51, 

CUN  T  I  NOE 

i  A' 

GU  TU  1929i 

bAl 

1  S2  6' 

I  FINAL  =  IFIl'lAL  ♦  1 

A  ^ 

)U  1  '>2  t>(  J  J  =  1  ,  N  S 

3A  > 

IKVIJJI  -  1  .  E  ♦  35  )  19280  ,19265, 1928C 

b  A  * 

IS  265 

I  f  (  At  (  l  ,  J  J  1  )  1  92  tie,  192  BG  ,1927'J 

b  AS 

1s27 

V(  JJ)  =  LL  (  I  ,  JJ)  -  Dt  (  1  ,  JJ)  *  um 

3  40 

i  S’c'. 

CUNTINUL 

*5  r 

1S2S' 

CliN  T  IN  Lb 

bbb 

GU  Tu  19220 

f. - UNT1MAL1JY  TEST  19300-19399 - 

3  b  3 

ISbU 

J  I F  M I  li  -  i  . 

3  bA 

3  u  1  9  j ,  2  1=1, NS 

»  k  I, 

IflUll)  -  KFMIM  193ol  ,  1  SJ,>2, 19302 

3bC 

1  S3 1 1 

UlFM  III  =  L(I) 

HI 

iL.JTtK  =  1 

it o 

JtNTEit  =  NSS 

bt  1 

i  s  j  .  ;• 

C  4. i  m  1  I  U  LL 

;•  1 3 

Jli  1  9  1 1  A  J  =  1  ,NS 

2  c  ^ 

IUVIJI  -  UIFMIM  19  3t  3,  1933a, 193  3  A 

btb 

1  S3'  3 

(J  1  F  rl  1  N  =  V(J) 

2  £  c 

J  t  IS  f  f:  K  =  J 

iti 

I  2  n  T t 1<  =  F  S  S 

3  ?v. 

1  S3  'A 

C  i  In  11  NLt 

3  7  2 

i)u  19)21  1=1, NS 

37b 

Ju  1 v b2  1  3=1  , N  S 

b  7b 

IMAM  l  ,  Jl  )  14  32  v,  ,193Gb,  19  32U 

37b 

1  S3.  5 

3  1  F  =  UL  (  1  ,  J  )  ♦  U(I)  ♦  V(J)  -  C  E  ( I  ,  J ) 

b  7 1 

IK  JIFI  19319,19320,19320 

3  7  7 

l  S3  1 1 

IFtuiF  -  oIFNIN)1931b ,19320,19320 

Ai  ( 

1S315 

UlFi-.  In  =  UIF 

A.l 

ILulhx  =  I 

A.  2 

Jr  MlLt  =  J 

A  1 

1  S3  21 

Ci.’ix  T  I  \ L> f. 

A  b. 

1  Sb21 

CUNTlNOt 

A'.  7 

IF  (  U  1 F  rl  I  iM  ♦  UELTA2)  19  Au  0,19325, 19325 
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FORTRAN  SOURCE  LIST  iAT MAX 

ISO  SIIOKCC  STATEMENT 

C - MI!  19>2i>  IS  FUR  OPTICAL  SOLUTIONS - ALL  CIF  ARE  NuN  NEGATIVE 

Alt  1  S 3 2 S  COOII  <OF 

All  L-.LL  I OU  AY  IltlTlHE.IUAT) 

412  Tli-it  =  FLOAT(1TIME)/60. 

A  1  A  Wti  Tl( c, ISl 15 (TIME 

A  1A  Rt'  1  I  E  (  o  1 1 9'  18)11 

Alb  ,Jhl  Tt  (  <  ,  19i  251UIFMIN 

a  1 1.  hi\  1  T  f.  ( t>  *  1  S*  .1 1  I 

a  1  7  RR  1  T.:  (  6 »  l  ‘U  1  7  ) 

A  2 Pr>  1  3AL  =  V. 

A21  90AL 

A  22  Jo  lUjA.  1  =  1, MS 

A  2  3  0  U  19  3  3  3  J  =  1 , N S 

A2A  IHALI  I  ,  J  >  II  9335,19331),  19330 

A  2  3  1  3  ;  31  TlMP  =  XX(  I  ,J)  *CF(  1  ,  J) 

A  2  t  P i  i  i*l A L  =  F'F  1  MAL  ♦  TEMP 

A 2  7  rtf  Hr.  (  o,  19,  2  1  )  I  ,  J,XX(  I  ,  Jl  ,  TFMP 

A'1',  1933 S  C  u  N  T  1  N  0  F 

A  2  <;  1S3AI  CONTInCC 

a3A  ftisl  It  I-i,  l*!..  19)  PRIMAL 

AiE  rt»l  IE  It.,  19  22) 

A  2c  00  1  >3l-  I  =1  ,MS 

A  3  7  WKl  Tf  (  t  ,  IS,  2  3)  I  ,U<  I  ) 

A  A  0  UUAL  =  OIJA  L  ♦  A  A  (  I  )  *  0(1) 

a  A  1  1  S  3  S  ■  LON T I  NCR 
A  A  3  Jo  lV3o.  J  = 1 ,  N  S 

AAA  R  R  I  I  E  (  o  ,  1 '»(  2  A  )  J  •  V  (  J  ) 

A  A  5  UUAL=  oOAl  +  (ill(J)  *  V  (  J } 

A  A  t  1  3  j  t  C  )i\  T  I  v  UE 
A  5  *»*  I  Tl  (  r ,  1 9.  2  ; )  DUAL 

a e i  call  exit 

C - KtPPE  SlnT/.  I  ION  OF  ENTERING  VECTOR  3Y  CURRENT  BASIS  1940o-19A99 - 

A3  2  1  SAL  i1  AE  AK  (  M  SS  I  =  O  . 

A  E  2  BdAn(jSS)=  i.  . 

A3A  F  S  P  A  C  c  =  i 

A  E  A  0(1  l  9 A  •  1  I  =1  ,MSS 

A  3  A  I  S  UiM  (  I  )  =  C 

A  7  AE  I  I  ,NSS  )  =  ANAK  (  I  ) 

At'  ISA.  1  3 1  (  I  ,MSS)=l  . 

Att  OF.  1  9*,l2  J-1,NSS 

At  2  J  SUti  I  J  )  =  l 

AtA  aL(MSS,J)  =  f)  D  A (  J ) 

A fj 3  ISa  2  Ut(MSSfJ)  = 

At?  JO  )9a.A  I=1,MSS 

A?  Ou  19A  A  J=1,NSS 

a  7  1  At'ASl  S(  I  ,JI  =  AF.  (  I  ,  J) 

A  ?  2  Y  Y  (  I  ,  J  )  -  v. 

a  7  1  IF  (  AE  (  1  ,  J  )  )  1 94i  A  » 1 9A0  A ,  19A  .'3 

A  7  A  ISA' 3  ISUM(I)  =  I  SUM  I  I  )  ♦  1 

A7E  JSUKIJ)  =  J  SUM ( J ) ♦  1 

Alt  N  SPACE  =  o  SPACE  ♦  1 

A7  7  1  SA 'A  CoNTlNUt 

5  w  1  1  SA..S  CONTI  Noe 

A.  i  I  1=1  1M 

U.  I  F  I  I  l  N  T  l  I<  -  MSS)194o7, 19416,19407 
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SOUkCiI  STATtMtNT 


F  Li  R  T  KAN  SOURCE  LIST  MATMAX 


I  S  M 

5.5 

1  S  4  (  t 

5  h 

5  .  7 

5  1  . 

1  54  .  7 

511 

1  54 1.8 

512 

5  13 
5  14 

1  54 1. 9 

51  = 

5  It 

1541. 

51  7 

15411 

5t 

5  £  1 

5  20 
53  J 
c  34 

15412 

542 

4)  ■* 

444 

15413 

552 

5  ;  j 

154  14 

5=4 

5  5  5 

5  5i 

15415 

5  5  / 

15417 

5  5  , 

K1 

15420 

5,-0 

5t  J 

1  5425 

5  r.  4 
=  c  *i 

1  54  3.. 

5  0  t 

1  5435 

7 

5  1  . 

571 

5  72 

1  544 < 

5  7 1 

5  74 

1  5445 

5  it 

1  545' 

5  77 

t  1 

t  1 

c  ,2 

1  54y5 

t  1.  4 

t  5 

1  4  4  0 

f  .  t 

0  7 

1  44  7' 

t  1  . 

Ml 

1  54  75 

t  12 

(  1  = 

1  4 

1  54C. 

t  1  0 

1  5  4  r  L 

1 1  7 

i  ' 

UF  (  I  L  .,T  L  •»  »F|  S  S  1  =i'  • 

)F(F>SS,J::N'TLh»  =  1. 

Gli  T I )  I'KI1. 

IMJrMlM  -  NSS)  IW)'),  19408,  19409 

U  F  (  ll.  llFr(,N.jSI=l» 

ul  MSS.JlUfr)  = 

liiJ  Tu  194li 

ill.  (  M  4 5 » JF  N To  K  )  =  I. 

Jt(  IcNTI:  <,nSS)  =  UE  I  ItNTfcR.JENTERI 

IH11  -  NluTAL)  19414,19414,19411 

(•KITtli-.lV)  ill)  IF  MIN 

WF1  Tc  ( (>  *  l  Gtj'j  3  I  I  1 

WhCl  U  (t>,  19'  161(1  SUM!  1  )  ,1  =1  ,MS1 

WK1  Tc(  6,  I9.ile>(  JSUMIJl  , J=1 ,NS) 

Du  17412  1=1,  MSS 

WRl  Ttlt.lS'i.  4)(At(  l  ,J>  ,J  =  1  ,NSS) 

Wrii  It  (6,  IV  l6)ItNTEK,JENTER 

Du  i  =  i »mss 

HKl  It  ( 15.  ',41  lAbASlbl  I  ,  J  )  ,  J=  1 ,  NSS  ) 

ul)  T  i  i  1 5 1 P  5 
Cunt inw 
M  I'JTflL  =  ' 

I  K  l  P  "•  r  •  =  i 

IFUTuIAl  -  MrUu>  19417, 19495,1949!) 
lKtt.\E  =  lKtP^fc  ♦  1 

IHUt-P-xC  -  y  U  I  G  )  1942,,  19420,19411 
UiJ  19455  1  =1  ,MS 

IHISun(l)  -  1)19455,19425,19455 

Ut)  1  O',', 3  J  =  1  ,NSS 

IF(  A  3 A S 1  !,(  I  ,J»  )  19445,19445,19430 
IMJ  -  NS  Silo  435,1  944 u  ,  19435 

Y  Y I  1  ,  .1 )  =  l  L  <  1  ,  N  S  S  )  /Ufc  (  I  ,  J  I 
JH  1  *  N  S  S  )  =  Y  Y  (  I  ,  Jl 
OL(MSS,JM  Ut(MSS,JJ  -  Y  Y  ( I  ,  J  ) 

111  lu  1 94  5| 

YY(  1  ,J)  =  t,L(  1  ,N  SSI 
1  j  U  TD  1  9  *,  t 
uUu  T  I  ,Ct 

I  SlMl  1  )  =  1  SIHC  I  )-  1 
JSU.-il  J)  =  J  S  U  M  (  J  J  —  1 
NID1AL  =  MuTAL  ♦  ! 

AoASl S( I , Jl=  0  » 

CimTI  mLE 

IFInIJTAL  -  Mbl  M  1946 ’J,  19*, 95, 19495 
DU  1  44", ■'  J=  1  ,  NS 

IK JSUMlJ)  -  1 1 19490, 19470, 1949U 

DJ  l  14  *  . •  I  =1  ,  MSS 

IF  UbASISl  I  ,  J)  119430, 19480, 19475 

Y  Y  (  I  ,  J  I  =  OF  MSS,  J) 

)L (  1 , NSS)  =  Ul ( I , NSS ) -  Ot(I,J>*  Y  Y  (  I  , J  ) 

vjD  T  I  19445 
C  i  j  4  T I  ,  CE 

I  s ui-i ( I )  =  isc"(i)-i 

J  SU'I  (  J  )  =  JSUM(J)-1 
Ail  A  S 1  SI  1  ,  J)  = 
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FORTRAN  SOURCE  LIST  MAT  MAX 


1  Si. 

(  2  1 
<  2l 
6  i  i 
/.  2  ~ 


tit: 
bl  I 
1 i : 
oji 
<  i  t 
MJ 
b  24 
(  j  4 
Kit. 
bit 
>  A 
b  A  1 
r  A  J 
l)i  i 
r.  A  A 
A  5 
A 


f>5  2 
i*  5  J 
r>  9  4 
t  ;  5 
(  56 
c  5  7 
o  6  1 
ft  2 
f.  f.  A 
fcf  i 
t  L 
f  7 
t  I  I 
f  7  2 
<•  7? 

6  7  5 
t.  7t 
o  7? 

7  ' 

7  l 
7  ; 
7  3 
7  A 
<<  j 
7 

7  ;  7 

7  1  . 
I  1  l 
7  1  /. 

7  I  i 


if  T Ll  IAL  =  NTUTAL  ♦  i 

1  449i  Cl/  f  I  I  HUt 

..I.  tl,  tjais 

1  54  9"  CtlN  I  1  fg  OF 

t. - VI.CliH  III  LtAVL  bASlS  IS  NUM  Ofc T t KP1  NfciJ - 195CC  -  17599 

155  .  THL  I A  =  1  .  h  f  j>5 

-)u  1,:>29  1=1, MSS 
UU  J  =  1  ,  rj  3  S 

IMAtl  l,J»  1195,7  ',  1  952",  195.4 
195  4  IF ( YY< l ,J) > 1 7520,1952', 195U6 
15526  !l  I  1  -  MSS)  1751',  19558, 1951( 

194  o  AL  TL  M.»  =  n  AAa  (  J  )  /  Y  Y  (  l  ,  J  ) 

,,U  T  Ij  l-)‘.lt 

l9bli  1FIJ  -  M5S)  1  7514 ,19512, 19514 
1.512  A  L  T  L  v  R  =  At- Ak(I) /YY( 1  ,J) 

Si )  k>  19  51  t 

19514  A  I.  T  t  M  R  =  X  X  (  I  ,  Jl/YYl  I  ,J) 

ISolc  IFULTMK  -  THLTA11951H, 1952  J,  1  752o 
1  S  5  1 ...  TiIlIA  -  /llr.MP 
iLcAvl  =  1 
JLlAVc  =  J 
15  521  Ct-JllfUh 
1  5  r>  2  5  C  ( ■  y  T  I  'vUl 

L - U-VfSF  J.-t*  j  l  Li  SLLUIICN  IN  T  t  R  7-3  Lf  CURMFINT  bASIS  1  tb  j  -  1  5699 

C - ,At»  Sol'JT  1  J\  IS  RLTURNFD  TC  192 jJ - 

lctf,  ju  lsj/-  1=1, MS 
Hj  1,.>1)  J  =  1  » N  S 
IF(Ar(  I ,JI  II  )  ,15 ,19615, 19605 
14-.S  IF  (  YYU  ,  J)  )  1  751-  ,  1  9o  1  5 , 1  76  lu 
l9tlL  XMl.J)  =  XXU,J)  -  THE  TA*YY  (  I  ,  J  ) 

14615  Cw,,  I  1  .JOl 
1  4 1"  2  J  C  u  4  1  1  f  Ct 

'  Ji  j  19')  1  •  1  =  1, MS 
IF  1  YY<  1  ,  sSS)  119o25  , 19o3 i"1 ,19625 
19625  A 1} 4 K  1  1  1  =  AtiA-i(l)  -  TrttTA  *  YYll.NSSI 
1  So  i  C  UN  T  I  HCo 

Jli  19*. a  /  J  =  l  , -M S 
IF ( YYlMSS, J)  )  19639 ,1964r ,1 9635 
IS'  I:  IMUkIJI  =  LHAK(J)  -  ImlTA  *  YY(MSS,J) 

1  4  o  A  i.  C . .  i ,  I  I  N  0 

AL  (  1  Lt.A  VL  ,  JLtAVt  )  =  r). 

UllLuAVi.  -  MSS )  1  969  2 ,19641,19  642 
19«)41  i)  u  A  R  (  JL  L  A  V[  )  =  J  . 

,ii  I  i  J  1  So  A  A 

191-42  !t(Jlt/-Vf  -  xSS  1  1  9644 ,1  464  3,  1  9644 

I  464  3  A.'ARl  ULAVtr)  = 

II  i:44  I  t  (  l  f  if  1 1  <  -  , MSS)  1  744  i  ,1  9644, 1965  ) 

19145  SHaoI-JI  i.lro)  =  TflfcTA 

a  f  (  s  s  s ,  j  i:  N  T  (  k  )  =  theta 
-,u  In  1 9  7i  ■ 

19  r.  5 '  l  f  1  J  L  Y  1 1  i.  —  NSS)1966',196A5,1  766‘) 

19o45  A,/M,v(  Il.-j  ff  i  )  =  THL  TA 

At  (  ItNTL^.isSil  =  THLTu 

G  'j  lit  1  '7  7 '  L 

lstot  xxi  i  :mh»  ,  jt \n.<i  =  imlta 
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FORTRAN  SOURCE  LIST  MATMAX 

ISN  SUUKCc  STATEMENT 

7  14  At(  ItNli.R,  JLNTER)  =  C  E  (  1  E  M  EK ,  Jt  NT  E  R  ) 

715  1KILLAVE  -  MSS )  l  96  75 , 1  97(  0 , 19675 

7  1 1  14675  IFIJLLAVt  -  NSS  >  19680  , 1  S7f  0 , 1 9  7l,0 
717  1  >6 60  XX(  iLtAVL .JLcAVE)  =  9. 

L - FrcgULNCY  Ilf  UbJ.  FUNCTION  PRINTOUT  CCNTKCLED  BY  SETTING  IrfRITE 

7  2.  1S7CC  [PRINT  =  1PKINT  ♦  1 

721  IFIIPkInT  -  l RKl TE) 19200,  197  10, 19710 

7*2  19  71U  I  PRINT  =  I 

723  DUAL  =  0. 

725  DU  1 9 7 5  i  1=1, MS 

72  5  1  S  7  6u  DUAL  =  DUAL*  AA(1)  ♦  U(  1  ) 

727  JIJ  197c.  J  =  1  ,NS 

730  1 S  7  f  '•  DUAL  =  DUAL  ♦  BB(J)  *  V(J) 

732  WR  l  T  u  (  6 , 1 91  39)11  ,  DUAL  ,1)  1 F  MI  N 

C - UtTAILtD  PRINTOUT  CGNTRCLEC  BY  SETTING  [TAB  =  1  IN  INPUT - 

73  3  1  9  60b  IMITAb  -  1)192  00,19810,19200 
734  1561C  RKl Tele,  1401  1) 

736  RR 1 T£(f ,19035) 

7  3  6  Wk I  If  I  6,19039) 

73  7  rml  TE  (  b,  19905)  (0(1), 1=1, MS) 

744  Rn  l  TE  (  6 , 1  9r  4u  ) 

74  5  RK  I  TC ( 6 , 19l05 ) ( V ( J) , J= l  ,  NS ) 

/  52  W.RI  TE(6,19i  3  3)  I  E  NTE  R  ,  JE  NT  E  R  ,  I  L  E  AV  E ,  JLEAVE,  THETA 

753  WRI TL(6, 19036) 

754  DU  1982^  1=1, MS 

75  5  WR1  TE(c,19(  27)  [ 

75o  l=b2i  Rkl Ttl6, 19105) (XX( I ,J» , J=l ,NS) 

764  MK1 Tt(6, 19037) 

7t5  RkI  Tt  ( t, ,  1  Ve  vd  )  (BBAR(  J)  ,  J  =  1  ,NS) 

772  WKI Tc( 6, 15C3B) 

7  73  RKl  T  E  (  6 , 1  9liC5  )  (  ABAR  ( I  )  ,1=1,  MS) 

ICOg  GU  TU  1921.0 

11'  I)  1  ENU 


Appendix  B 


SAMPLE  PROBLEM 


EXAMPLE 

A  short  example  is  given  here  to  illustrate  (a)  the  type  of  problem  to 
which  the  algorithm  can  be  applied  and  (b)  the  details  of  the  solution  process. 
The  example  with  some  modifications  is  taken  from  Ref  3.  The  problem  is 
one  of  allocating  several  types  of  commercial  aircraft  to  various  routes  (e.g., 
New  York  to  Dallas)  in  order  to  maximize  overall  profit  (revenue  less  opera¬ 
tional  costs)  for  the  system  (see  Fig.  Bl).  The  problem  can  be  stated  as 
follows: 

Let  m  denote  the  total  number  of  routes 

n  denote  the  total  number  of  types  of  aircraft 

a,  denote  the  anticipated  number  of  passengers  on  route  i  per  month 
b(  denote  the  number  of  aircraft  of  type  j 

denote  the  monthly  passenger- carrying  capacity  of  aircraft  type  j 
on  route  i 

r,  denote  the  revenue  per  passenger  on  route  i 

denote  the  monthly  operational  cost  for  operating  aircraft  type  j 
cn  route  i. 

Then  we  seek  to  find  the  allocations  i|(  of  aircraft  type  j  to  route  i  for 
the  system  that  maximize  profit  subject  to  constraints  on  the  number  of  antici¬ 
pated  passengers  using  the  various  routes  and  the  available  number  of  each 
type  of  aircraft. 

Specifically 
maximize  profit  r(x(j ) 
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Route  i 

Aircraft  type  j 

1 

2 

3 

4 

i 

1600 

900 

2 

1500 

1000 

500 

1100 

3 

2800 

1400 

2200 

4 

2300 

1500 

700 

1700 

5 

8100 

5700 

2900 

5500 

o.  Aircraft  Passenger-Carrying  Capocity  par  Month  (d-) 


Route  i 

Aircraft  type  / 

1 

2 

3 

4 

i 

18,000 

17,000 

2 

21,000 

15,000 

10,000 

16,000 

3 

18,000 

16,000 

17,000 

4 

16,000 

14,000 

9,000 

15,000 

5 

10,000 

9,000 

6,000 

1,000 

b.  Oparationol  Costs  in  Dalian  par  Month  (s)(J 


Fig.  Bl— Sample  Data  for  Aircraft  Allocation  Problem 


Blanks  in  cell  /,  j  indicate  that  aircraft  type  /  is  never  assigned  to  route  i. 
*»  5  (routes)  n  =  4  (typos  of  aircraft) 

Passenger  data  Passenger  fare  doto  Aircraft  data 


=  25,000 

r,  =  $130.00 

o- 

II 

o 

12,000 

r  2  =  $130.00 

b  J  =  19 

=  18,000 

r3  =  $  70.00 

f>3  25 

9,000 

r4  =  $  70.00 

b4  =  15 

=  60,000 

r5  =  $  10.00 

n  i) 


operational  costs 


£  1  (r  d  -  s  )x 

,  1  I  ,  •  *J  »l  -J 


subject  to 


i.  d  x  la 
,  ,  «M/  » 


(l  I . ffl)  do  not  exceed  passenger 

demand  for  each  route  I 


,V-b- 


ri)  do  not  exceed  aircraft 

availability  for  each  type  j 


(Bl) 
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Data  for  m  =  5,  n  =  4  are  given  in  Table  1  as  taken  from  Ferguson-Dantzig.3 
It  should  be  noted  that  if  c(.  =  -  s1(  above,  then  the  form  of  the  problem 

is  that  of  Prob  3.  The  Ferguson-Dantzig3  example  requires  that  all  aircraft 

m 

be  allocated,  i.e.,  .1  =  b(  (j  =  1 . n).  The  algorithm  here  does  not  require 

that  all  resources  be  allocated.  However,  at  least  all  the  row  or  all  the  column 
constraints  will  be  binding  for  an  optimal  solution.  If  all  the  column  constraints 
are  binding,  then  of  course  the  equality  constraints  on  the  columns  are  satisfied. 

The  example  of  Ferguson-Dantzig3  can  be  formulated  from  the  example 
here  for  m  =  5,  n  =  4  as  follows: 

Let 


then  we  have 


subject  to 


*.s 

s  1 5  =  ri  d  i5  "  1 


(B2) 


max  zCl, .)  =  max 

“‘I  '<) 


5  5/4 

S  r .(“  -  i  .)  -11  s.jt. 
i=l  1  1  15  .=  1  (=1  "  ‘I 


max 

X- 


5 


1  r  a  -  min 
'=1  '  x 


5/4 

1  1  s  (x  +  r  x  el 

t=l  )=1  ,J  '>  1  ,5I 


5  /  5 


‘-Mi-1 


(B3) 


J.V.I 

(i  =  1,  . 

.  .  4) 

4 

,=Y'i =  b< 

()=!.•■ 

.  ,  5) 

Ferguson  and  Dantzig3  consider  the  problem 


5,5 

*1 

1=1 


subject  to  the  two  sets  of  constraints  given. 

Details  of  the  solution  are  printed  at  each  step  of  the  solution  for  the 
problem  expressed  by  Eq  Bl.  A  completely  detailed  description  illustrating 
the  algorithm  on  a  step-by-step  basis  is  presented  between  the  first  and  second 
intermediate  printouts  (iterations)  of  phase  2  of  the  algorithm. 
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INPUT  CONSTANTS 


INPUT  CONSTANTS  FCH  OPTIMAL  ALLOCATION  PRUULtM 


MS=  5  NS=  4 


INPUT  VALUfcS  FUk  CF (  I  ,  JI 


RUW 

1 

L. 19000 JE 

06 

0. 

C  . 

0. lOOOUuE 

06 

ROW 

2 

_  0. 1740cOE 

Oo 

o.  1  ISC  CUE 

06 

0.550000E 

05 

0. 127000E 

06 

ROW 

3 

L.  1  6  BOOut 

06 

C.  B20C  Ou  E 

05 

0. 

0.1370U0E 

06 

ROW 

4 

2.14:0006 

06 

u . 9 luOOC  E 

05 

0.4O00UCE 

05 

0.104000E 

06 

ROW 

5 

_  _  0. 710UUOE 

05 

o .  4BO(!OOE 

C  5 

0. 230)006 

05 

0.450000E 

0  5 

INPUT  .VALUfcS 

FUR 

UE(  I  ,J) 

ROW 

I 

_____  C. 16CJCJE 

04 

C.  100000E 

01 

0. 1C 00 00 E 

_01 

0 . 90000J  E 

03 

ROW 

2 

C. 15O00UE 

04 

0.  l'JOCOOE 

04 

0. 500000 E 

03 

0.1 10OU0E 

04 

ROW 

3 

0. 280000E 

C4 

0. 140000E 

04 

0. lOOOCOE 

01 

0. 220000  E 

04 

ROW 

4 

0 . 2  3000 JE 

O 

0. 150t  JOE 

04 

0. 700000 E 

0  3  _ 

0 . 1_7QOOO  E 

04 

ROW 

5 

C.aiOJ'jJE 

04 

C . S7COCOE 

04 

0. 29000C  E 

04 

0 • 5  50000  E 

04 

INPUT  VALUES 

FOR 

A A ( I )  OF  THE 

RCW 

CCNSTRAI NTS 

0.25UUG Jt 

06 

C.  1200C0E 

05 

0. 180000  E 

05 

0. 900000E 

04 

INPOT  VALUES  FUR  tJH  C  J 1  OF  THE  CULUMN  CONSTRAINTS 


C.lOJQoJt  u 2  0.19(0000  02  0.250000E  02  0.15U00oE  02 


FREO.  UF  OBJ.  FUNCT.  PRINTOUT  1 


1--1  4.  4^-4. 


:h  OPTIMAL  ALLOCATIUN  PRDULfcM 

tfcu.  (JF  OBJ.  FUNCT.  PRINTOUT  1  DETAILED  PRINTOUT  IF  I  TAB 


C.  0. 100000 E  06 

0. 550000E  05 _ 0.127000E  06 _ 

0.  0. 13  70UOE  C 6 

0  •  41/0000  E  05  0.104000b  06 

0. 2  30  H  OE  05  0.450000b  05 _ _ _ 

G.lCJOOOEJJl  0.90Q00JE  03 _ 

0. 500000  E  03  0.110000E  04 

0.1 00000  E  01  0 . 220000  E  04 

0. 700000 E  03. _  0.1 70000 E  04 _ 

0. 29000C  E  04  O.5500OOE  04 

CONSTRAINTS 

0 • 180000  E  05  o7sOOOOOE_C4  0.600000E  05 
MN  CONSTRAINTS 

0.1 5U00o  E  02 _ 


0.250000E  02 


PHASE  1,  ITERATIONS 


MATKIX  MAXIMUM  ITcKAIluN  NUMBER  l 

ALLOCATION  SELeCTEJ  rJW  I  CUL.  I  XXII, J)=  G.10000E  02 

VALUE  UP  OBJECTIVE  PUNCTI UN  =  O.19G00E  07 

MATKIX  MAXIMUM  1  It  K  AT  I  IN  NUMBER  2 

ALLUCATIJN  SELECTED  RUM  3  CUL.  4  XXII, J)=  0.818I8t  Ul 

VALUE  OP  OBJECTIVE  FUNCIIUN  =  0.302U9E  07 


MATRIX  MAXIMUM  ITERATION  NUMrfEK 
ALLOC A  I  I  UN  SELtCTED  KUw 
VALUl  Ul-  UBJECT1VE  FUNCTION 


3 

2  CUL.  4 
0.38b68E  07 


XXII, J)=  0.68I82E  01 


MATRIX  MAXIMUM  ITEkATIUN  NUMBER 
ALLOCATION  StLECTfcU  RUW 
VALUE  OF  OBJECTIVE  FUNCTION 


4 

2  CUL.  2 
0.44043E  07 


XXII, J)=  0.45000E  01 


MATRIX  MAXIMUM  ITEFtATIUN  NUMbFK  5 

ALLOCATION  SELECTED  KJH  4  COL.  2  XXII, J 0.60000E  01 

VALUE  OF  OBJECTIVE  FUNCTIuN  =  0.495G3E  07 

MATRIX  MAXIMUM  ITEKATIUN  NUMBER  b 

ALLUCATIJN  SELECTEO  kUk  B  CUL.  2  XXII, J)=  0.85000E  01 

VALUE  OF  OBJECTIVE  FUNCTIUN  =  0.535B3E  07 


MATRIX  MAXIMUM  ITEKATIUN  NUMBER  7 

ALLOCATION  StLECTEU  RUW  5  COL.  3  XXII, J)=  0.39828E  01 

VALUE  OF  UBJECTIVE  FUNCTIUN  =  0.S4499E  07 

ITERATION  1  PRIOR  VALUE  OF  OBJ  FUNCTION  0.54499E  07  MAX.  VIOLATION  0 


Jl=  0 • 10000 E  02 


J)=  0.818l8t  01 


J)=  0.68182E  01 


J)=  0.450O0E  01 


J)=  0.60000E  01 


J>  =  0.85000E  01 


J)  =  0.39828E  01 

54499E  07  MAX.  VIOLATION  OF  DUAL  CQNSTR.  -0.96428E  05 


UET4ILtD  I  NT  fcRMEDI  AT  E  PRINTOUT 


PklUK  VALUES  UF  THfc  DUAL  VARIABLES  U(I)  IN  OROfcR 

0.  0.1122O7E  03  0 .606489 E  02  0.588046E  02 


PRIUR  VALUES  UF  THE  OUAL  VARIABLES  VIJ)  IN  ORDER 

G.lSuOOOE  06  U.27V310E  04  0.  0.357242E  04 

ENTERING  VECTOR  I  4  LEAVING  VECTOR  2  4  XXIIfcNTER.J 


CURRENT  SOLUTION  ARRAY  XXI  I, J) 
ROW  1 


ROW  2 
ROW  3 
ROW  4 
ROW  5 


C . ICO  300b  02 

U. 

c. 

c  • 

0. 

UNUSED  RESOURCES 

C. 

UNUSED  CAPACI TIES 


U. 

C .  1200E0E  02 
l)  • 

0. 6DU000  E  01 
O.IOOCOOE  31 


C. 

C. 

0. 

0. 


0.6818 18 t  Cl 


0. 

0.818  182  E  01 

0. 


C.187241E  02  0. 


COLUMNS  I  THRU  NS  IN  CRDER 
0.  0.627586E  01 

KUWS  1  THRU  MS  IN  CRDER 


0. 


ITERATION 


C.28S364E  .4  C.  0.  0. 

2  PRIOR  VALUE  OF  OBJ  FUNCTION  0.61074E  07 


MAX.  i 


1 1  L  t  L)  I  NT  ERMEOI  AT  E  PRINTOUT 


IAL  VARIABLES  U(I»  IN  ORDER 

22C7E  )3  0 .606489 E  02  0.568046E  02  0.793103E  01 _ 

AL  VARIABLES  V(J)  IN  ORDER 

931CIE  04  0.  0 . 35  7242  E  04 

EAVING  VECTOR  2  4  XX ( I  ENT EK , JENTER  I  =  0.68182E  01 


> 


C. 

OOLOE  02  G. 

0. 

0000 E  01  0. 

DC DOE  31  0.187241E  02 

UMNS  1  THRU  NS  IN  CRDER 

0.6275 86 E  01 
KS  1  THRU  MS  IN  CROER 
0. 

DF  OBJ  FUNCTION  0.61074E  07 


0.681818 E  Cl 
0 .  " 
0.818 182  E  Cl 
0. 

0. 

0. 


0. 


MAX  . 


0. 

VIOLATION  UF  DUAL  CONSTR. 


-0.55661E  05 
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DETAILED  DESCRIPTION  OF  ONE  ITERATION  OF  PHASE  2 


A  detailed  numerical  examination  of  one  iteration  of  phase  2  of  the  algo¬ 
rithm  is  presented  here.  Given  the  solution  array  XX(I,J)  of  the  previous  page, 
we  wish  to  test  optimality  of  the  solution.  Proceeding  as  in  step  2  of  the  algo¬ 
rithm  we  solve  Eqs  15  for  (u, . u-.  i»| . r,  i. 

Since  xr)}  =  6.27586  >  0  and  xr)t  =  2863.63  >  0,  we  have  i  t  =  0  and  ut  0 
immediately,  i.e.,  we  solve  for  the  zero- valued  dual-space  variables  first. 

The  remaining  seven  equations 


< i 


m  •  r 


i  i 


if 


X 


I) 


0 


are  solved  sequentially.  We  have 

xn  10  h=i  i  fii 

t ,  ,  6.81818  0  =  r  ,  , 


1 90.000 
100.000 


2:1.000 

2.000 


7.0:110.1 


1  2 

3  ., 
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1  ()  =  !■_,  11.  tH.ooo  i:,;oo)(:.o:iiui)  2701.10 


0  =  11., 


».,2  >' 


0  =  u 


1*2  1  2 

*  12  1  2 


7.7)  l.'il."  0  =  11, 


if  4  2 

C  U  1  1 


‘31 


I  I ',.000  270:1.10 
1000 

91.000  270.1.10 

17,00 

1 17. (KM)  100.000 

2200 


112.207 


,8.80  If, 


16.8182 


Moving  on  to  step  3  of  the  algorithm  the  optimality  test  is  now  made  (i.e., 
Conditions  7D  in  scalar  form  of  the  m  +  n  +  rroi  =  29  equalities  of  Eq  16  are 
tested).  It  is  seen  directly  that  u,  («  =  1,  ....  5)  and  r  (j  =  1,  . . . ,  4)  are 
nonnegative.  Seven  of  the  remaining  mu  =  20  inequalities;  namely,  «s  = 

Vi  +  ~  f.,  =  0  for  'ii  >  0  (i-e.,  for  x,„  x,„  *S2.*22 .  and  *:u>  are 

equalities,  hence  only  the  remaining  20-7  =13  values  of  us  are  tested  for 

'  0 

nonnegativity.  Computing  directly  we  have 


ir 

"\2 

dl2“l 

*  I’o 

1  |2  (IHOI  .  279.1.10  0 

0 

w 

•'ll 

1  3I(  1 

‘  1  3 

c  n  (Dio)  .  0  0 

0 

14* 

s2 1 

J2|«2 

'  '  1 

1  ,,  (17,0(0(1  12.207)  •  190.000  171.000 

0 

ir 

s2.1 

3 

f  *3 

i  21  (7,0(0(112.207!  .  0  7,7,.()(K) 

0 

35 


RAC 


ui  *  d24u2  +  "4  -  <24  -  (1100)012.207)  +  100,000  -  127,000  >  0 

u;  »  d3,u3  +  i>,  -  c3I  =  (2800)06.8182)  +  190,000  -  168,000  >  0 

ws  -  d32u3  *  v2  ~  c 32  *  (1400)06.8182)  +  2793.10  -  82,000 

32  =  -55,661.42  <  0 

IB  =  ^33m3  f  ^3  -  c33  =  0)06.8182)  +  0  -  0  >0 

«;  =  d4|u4  +  u,  -  c4,  =  (2300K58.8046)  +  190,000  -  145,000  >  0 

W5  =  d43u4  +  i>3  -  c 43  =  (700X58.8046)  +  0  -  40,000 

43  =  41,163.22  -  40,000  >  0 

ki  =  d44u4  +  ^4  -  c44  =  (1700X58.8046)  +  100,000  -  104,000  >  0 


w  *  d5jM5  ♦  i>,  -  f5,  =  (8100X7.93103)  +  190,000  -  71,000  >  0 


ws  =  d5iu5  *  "4  -  <54  =  (5500X7.93103)  +  100,000  -  45,000  >  0 

54 

It  is  seen  at  this  point  that  the  only  violation  and  hence  the  maximum  vio¬ 
lation  (step  4  of  the  algorithm)  of  Conditions  7D  of  Eq  16  is 


w  =  -55661.42  <  0 
s32 

Hence  the  solution  is  not  optimal  and  the  vector  d32?3  +  e7  enters  the  basis. 

The  vector  selected  to  leave  the  basis  is  determined  next,  but  first  it  is 
necessary  to  determine  the  components  of  the  entering  vector  d32?3  +  ?7  in 
terms  of  the  current  basis.  The  set  of  m  +  n  =  5  +  4  =  9  linear  equations 
Eq  21  is  thus  solved  (step  5  of  the  algorithm).  We  have  for  the  case  at  hand 
the  following  system: 

y  1 1  (1600)  +  yf* (900)  +  y 32 (1)  -  0 

y|f(1000)  -  0 

y\l (2200)  -  d32  -  1400 

yj| (1500)  =  0 

y^  (5700)  +  y32  (2900)  =  0 

v32  -  0 
>11  -  u 

>51  *  >21  ♦  >22  -  > 

>53  *  >63  -  0 

’ll  ♦  >55  -  o 
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The  solution  to  the  foregoing  equations  is 


32  _  o  32  .  HOO  32  _  0  32 
>22  -  u-  y34  -  2200'  >42  U>  > 


52 


>•  >n 


o, 


v32 

>53 


5700  .32 

2900  f>63 


5700  32 

2900  '  y 


1400  32 

2200’  y'5 


900)(  1400) 
2200 


Thus  the  entering  vector  d32?3  +  ?7  =  1400c3  +  e7  can  be  written  as  the  follow¬ 
ing  linear  combination  of  the  current  basis  vectors. 

<  900M  1  400), 


1400?,  .  ?, 


2200  |90L^I  4  *  2200 

+  ^  (2200?,  ♦  ?9)  ♦  115700?,  ,  ?.l 


-[?,] 


5  l  00  ;  — ♦  I  5  1  IfU  1 

2900  2  ^5  *  c8'  "  2900  ^8 


5700 


The  vector  now  selected  to  leave  the  basis  is  determined  from  the  indexes 
that  yield  the  minimum  in  the  brackets.  Using  Eq  29  (step  6  of  the  algorithm) 
we  have 


8.18182 

1  6.27586 

2863.64 

min 

1400 

'  1  ’  /  5700\ 

’  (900X1400) 

2200 

\2900; 

2200 

=  l 


and  the  minimizing  indexes  are  (5,2),  since  -Is  =  yields  the  minimum  value 

yS2 

0 .  Thus  the  vector  d52?5  +  ?7  =  5700e5  +  e7  leaves  the  basis. 

The  new  solution  is  evaluated  (step  7)  using  Eqs  28  and  appears  in  the 
next  intermediate  printout.  The  algorithm  again  returns  to  step  2  for  the  next 
iteration.  The  value  of  the  objective  function  was  increased  by  §i(-u’  )  = 

$55,661.42  during  this  iteration. 
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MC 


PHASE  2,  REMAINING  ITERATIONS 


DE  TAILED  I  hi  ERMEDI  AT  b  PRINTOUT 


PrIOR  VALUcS  uF  1 Hfc  DUAL  VARIABLES  U(I>  IN  ORDER 

H.1I22U7E  (3  0.168182E  02  0. 588^461 

PkIUR  VALUES  Uf  iHt  DUAL  VARIABLES  V(J)  IN  ORDfcK 

C.lVOK'it  Do  l'  .2793I0E  J4  U.  J.100UU0I 

ENTERING  VECTOR  3  ^  LEAVING  VECTOR  b  2  XXIIEl 


CURRENT  SOLUTION  ARRAY  XX!  I, J) 


KUm 

1 

C.u;  / v C E  J2 

0. 

0  • 

U . 74545b j 

ROW 

2 

L  . 

f.l20GClE  j2 

C. 

‘  0.  * 

ROW 

3 

V.  • 

o.ICOOUoE  01 

0. 

0 . 7S454bl 

ROW 

4 

c. 

C  .60(3000  E  01 

u . 

0. 

ROW 

5 

G  • 

0. 

0.206897E  02 

0. 

UNUSED  RESOURCES 

COLUMNS  1  THRU 

NS  IN  CRCER 

C. 

C  . 

C.431034E  01 

0. 

UNUSED  CAPACITIES 

ROWS  1  THkU  MS 

IN  CROER 

- 

O.^VOVlt  04 

L  . 

0. 

u. 

ITERATION 

3  PRlUR 

VALUE  CF  OBJ  FUNCTION  3.61630E 

07  M 

Bailed  intermediate  printout 


)UAL  VAKlAbLFS  U ( I )  IN  URDER 

ll22u7F  I  0.168182E  02  0.588u4t>E  C2  O.  793103E  01 _ 

»U7L  VARIABLES  VIJI  IN  ORDER 
{ 7 3 1 0 E  14  0.  J.IOC'UUDE  06 

LEAVING  VECTOR  b  2  X X ( I E NT Ek , J EN T EH  I  =  O.IOOOOE  01 


J) 

o  • 

20CEIE  j2  C. 

X’OUoE  J1  0. 

)0000E  01  0. 

C.2G6897E  02 
.UMNS  1  THRU  NS  IN  CRCER 
C.431034E  01 
ms  1  THRU  MS  IN  CROER 
0. 

CF  OBJ  FUNCTION  D.61630E 


U .  74545b E  C 1 

'  0. 

0 . 7b454bE  Ul 

0. 

0. 

0. 

U.  0. 

07  MAX.  VIOLATION  OF 


OUAl  CUNSTR.  -0.26727E  05 


PR  lUR  VALUE  S  UE 

0. 

DETAILED  INTERMEDIATE  PRINTOUT 

THE  DUAL  VARIABLES  U(l)  IN  ORDER 
0.565455E  02  0.168I82E  02 

0.21 6970 E 

-i 

C  2  0.7931 

PKIUR  VALUES  OF 

THE  DUAL  VARIABLES  V(J»  IN  ORDER 

- 

0.1900UUE  06 

0. 584545E  05 

0. 

0. 100000E 

C  6 

ENTERING 

VECTOR  2  3 

LEAVING  VECTOR  6  3 

XXI IENTER, JENTER  I  ^ 

_  _ _ _ j 

CURRtNT  SOLUTION  ARRAY 

XXC I , J» 

T) 

ROW 

1 

— 

C. ICOOcOE  02 

0. 

0. 

0.8826U2E 

01 

ROW 

2 

0. 

li.  984483E  01 

0 . 4  3 10  34  E  01 

0. 

ROW 

3 

C. 

0.315517E  01 

r. 

0.617  398  E 

01 

ROW 

A 

C. 

G.600000E  01 

0. 

0. 

ROW 

5 

c. 

0. 

0.206897E  02 

0. 

— 

LNUSEO  RESOURCES 

columns  i  thru  ns  in  order 

— 1 

C. 

0. 

0. 

0. 

UNUSED  CAPACITIES 

ROWS  1  THRU 

MS  IN  ORDER 

— - -  — 1 

— J 

U.105658E  04 

0. 

0. 

0. 

0.  J 

ITERATION 

4  PRIOR 

VALUE  OF  OBJ  FUNCTION  0.62782E 

07  MAX.  VIOLATIOl 
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PRINTOUT 


IN  GROER 

B2E  02  0.216970E  C2  0.793I03E  01 _ 

IN  ORDER 

0.100000E  C  6 

3  XX( IENTEK«J ENTER  I  =  0.431C3E  01 


0 .882602  E  01 

3^e  oi  ~~o: 

0.617398E  01 

h  °* 

97E~02  '"O'. 

(DRCER 


EP 

0.  0. 

62782E  07  MAX.  VIOLATION  OF  DUAL  CONSTR.  -0.12853E  01 


DETAILED  1 NTEPMEDl ATE  PRINTOUT 


PRIOR  VALUES  OP  THE  DUAL  VARIABLES  U(I)  IN  ORDER 

0.  0.565455E  02  G.168182E  02  0.2169 


PRIOR  VALUES  OF  THE  DUAL  VARIABLES  V(J)  IN  ORDER 

L.19000UE  06  0. 584545E  05  0.267273E  05  0.1000(1 

ENTERING  VECTOR  5  5  LEAVING  VECTOR  l  5  X X ( 


C . 1 COOOOE  Ut 


CURRENT  SOLUTION  ARRAY  XX(I,J) 
ROM  1 

ROW_  2 _ 

ROM  3 
ROM  4 
ROM  5 


0. 

0. 

0. 

c7 


o. 

G.8G0000E  "6T 
0.500000E  01 
0 . 600000E  01 

o.~" 


0. 


O.iOOOOl 


UNUSED  RESOURCES 

c. 


UNUSED  CAPACITIES 

0. 


0.800000E  01 

C. 

0. 

C. 170000E~02 
COLUMNS  I  THRU  NS  IN  CRCER 
0.  0. 

RUMS  1  THRuTmS  IN  ORDER 

U.  0. 


0. 

0.500001 


0. 

or 


0. 


TOTAL  EXECUTION  TIME  FOR  ALGORITHM  =  0.7333333  SEC. 

NUMBER  OF  ITERATIONS  AFTER  INITIAL  SOLUTION  4 

MAXIMUM  VIOLATION  OF  DUAL  CONSTRAINTS  0. 
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Il  Lfc  D  1  NT  E  Ri-'EDI  AT  E  PRINTOUT 


AL  VARIABLES  U ( I  I  IN  ORDER 

5455E  02  0.168182E  02  0.216V7aE  02  -U.128527E  01 _ 

AL  VARIABLES  V  I  J  >  IN  ORDER 

4545E  05  0.267273E  05  U.IOOOUOE  06 

EAVINO  VECTOR  l  5  XX ( I  ENT ER » JENT ER  ’  =  0.10700E  05 


;  0.  O.IOOOODE  02 

)oooe  "ol  XT.  aoodoom  ol 

)000E  01  C.  0 . 500000 E  01 

)000E  01  0.  0. 

C . 1 70  0  DOE ”02  0 . 

JMNS  1  THRU  NS  IN  CRCER 

0.  0. 

IS  1  THRU"  MS  IN  ORDER 

0.  0.  0.107000E  05 

0.73333  33 _ S  EC. _ 

ION  4 

0.  _ _ 
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DtTAILS  CF  THE  SOLUTION 


ROW 

I 

CUL  UMni 

1  ALLOCATION 

0.  no  )0t 

02 

RETURN 

FRUM 

ALLOCATION 

( 

KUW 

L 

CLILUMN 

A  ALLOCATION 

0.  10000  E 

j  2 

RETURN 

FROM 

ALLUCAT ION 

I 

HUH 

2 

COLUMN 

2  ALLOCATION 

C .  8000GE 

01 

RETURN 

FROM 

ALLOCATION 

ROW 

2 

Column 

3  ALLOCATION 

O.aoOuOE 

0  1 

RETURN 

from 

ALLOCATION 

I 

KOM 

3 

column 

2  ALLOCATION 

0 • 500 00 E 

J1 

RETURN 

FROM 

ALLUCAT  ION 

< 

RUM 

3 

COLUMN 

A  ALLUCATIUN 

C.50O0UE 

01 

RETURN 

FROM 

ALLOCATION 

( 

KUW 

A 

COLUMN 

2  ALLOCATION 

O.6O0OOE 

01 

RETURN 

FRUM 

ALLOCATION 

ROW 

5 

COLUMN 

J  ALLOCATION 

C . 17U00E 

02 

RETURN 

FRUM 

ALLUCAT ION 

VALUE 

OF  PRIMAL 

LSJeCTIVE  FUNCTION  c  .6  29  200  E  07 

- 

VALUES 

OF  THE 

CUAL 

SPACE  VARIABLES  1 LAGkANGE 

MULTIPLIERS  , 

SHADOW  PRICES) 

- 

KUW 

1 

U<  I  )  = 

O.lJ.Itt  02 

ROW 

2 

U<  I  )  = 

u*64'vCw'E  l  2 

KUW 

3 

um  = 

U.221A3E  02 

RUW 

A 

U(  I  )  = 

C.26667E  02 

4 

RUW 

5 

UI  I  »  = 

0. 

COLUMN 

1 

VI  J)  = 

0.16917E  06 

m 

COLUMN 

2 

V<  J »  = 

j.blOCCt  05 

j 

COLUMN 

3 

VI  J )  = 

0.23UU-E  05 

COLUMN 

A 

VI  J)  = 

u.tt82F6E  05 

-i 

■J 


H 

-I 


VALUE  OF  CUAL  UEJECTIVE  FUNCTION 


0.629200E  07 


N 

FROM 

ALLUC AT IUN 

0.  190C0E 

07 

N 

FROM 

ALLUCAT ION 

0.  lOOCOE 

07 

N 

FKUM 

ALLOCATION 

C • 92CC0E 

06 

H 

FROM 

ALLOCAT 1UN 

0.44000E 

06 

N 

FROM 

ALLUCAT ION 

O.410C0E 

06 

N 

FROM 

ALLUCAT ION' 

0 . 68500E 

06 

N 

FROM 

ALLOCATION 

0.54600E 

06 

a 

FROM 

ALLUCAT ION 

0.39100E 

06 

(RiCkS) 


Appendix  C 


COMPARATIVE  SIMPLEX  ALGORITHM 

This  appendix  contains  a  general  simplex  algorithm  programmed  spe¬ 
cifically  to  solve  Prob  3.  The  program  accepts  the  same  input  data  as  the 
MATMAX  algorithm  with  the  exception  of  the  IWRITE  and  ITAB  information 
used  in  the  MATMAX  algorithm.  The  algorithm  here  has  been  used  for  com¬ 
parative  purposes,  with  particular  attention  given  to  the  times  required  by  the 
two  algorithms  to  solve  the  same  problem. 
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FORTRAN  SOURCE  L  1ST 

ISN  SOURCE  STATEMENT 

0  $  IB  FTC  MATLIN 

1  DIMENSION  CE(20,24) ,DE(2I,251 , AA ( 20  I , BB (24 1 

2  DIMENSION  AS ( 42  ,474 )  ,BS( 45  1  ,CS I  52 5 1 

3  DIMENSION  1  PA THS < 45 ) ,C TS ( 45 ) , BTS (45 ) , ATS I!45 ) 

4  DIMENSION  ZS(525) ,ZMCS<525 1 ,ATS J 1 52 51 

5  18G01  FORMAT ( IH1 I 

6  18002  FORMAT ( 1H  I 

7  18003  FURMAT (3110) 

10  18004  FORMAT ( 6F 12 • 6 ) 

11  18005  FORMAT (15X,6E16,6) 

12  18006  FURMAT(bX,2HMS,8X,2HNS,2X,15HPRlNT  FREUUENCY) 

13  18007  F0RMAT(30X,46H INPUT  CONSTANTS  FOR  OPTIMAL  ALLOCATION  PROBLEM//) 

14  18008  FURMAT( 15X»24H INPUT  VALUES  FOR  CE(l,J)//> 

15  16009  FORMAT! 15X.24HINPUT  VALUES  FOR  DE(I,J)//> 

16  1  do  10  FOKMAT( 15X,45H1NPUT  VALUES  FOR  AA ( I )  OF  THE  KUW  CONSTRAINT S/ / ) 

17  18011  FORMA  T ( 1 5X ,48HI NPUT  VALUES  FOR  BB(J)  OF  THE  COLUMN  CONSTRAINTS//) 

20  18015  FURMAT(//37HTOTAL  EXECUTION  TIME  FCR  ALGORITHM  =  , F 12 . 7, 1 X, 4HSEC . ) 

21  18019  FORMAT! //34HVALUE  OF  PRIMAL  OBJECTIVE  FUNCT IUN, E2G .6 1 

22  18027  FORMA T ( 6 X « 4H  ROW, 13) 

23  18020  FORMAT! 101 10) 

24  18028  FORMAT ( 6X , 7H  COLUMN, 13) 

25  MRI TE ( 6 , 1 8001 t 

2fc  l.'Rl  TE!  6, 18007  ) 

27  READ! 5,18003) MM,NN 

32  MRITE!6, 18006) 

33  MR 1TEI6, 16003) MM ,NN 

34  WRITE(6, 13002) 

35  WRITE(6, 18008) 

36  DO  18035  1*1, MM 

37  RE AO (5, 18004) (CE( I , J) » J*1 ,NN) 

44  MR1 TE ( 6, 1 802  7 )  I 

45  18035  WRITE(6,18o05) (CE( 1  ,  J )  «  J*1 «NM 

53  MRITE(6,1 800?  ) 

54  WRITE(6, 18009) 

55  DO  18040  1=1 , MM 

56  RE AD (5,18004) ( DE ( I , J I , J  =  1  ,  NN) 

63  MRI TE! 6, 16C27) I 

64  18040  MRI TEI6, 180051 (0E( I ,J) ,J  =  1  ,NN) 

72  MkITE(6,180o2) 

73  MR1TE(6,18C10) 

74  READ(5,18C04) < AA! I ) ,1 =1 ,MM) 

101  MRI  TE!  6, l 8005) (AA! II ,1 =1 ,MM) 

1G6  MRITEI6, 18002) 

1C7  MRITEI6, 18011) 

llu  READ! 5, I80C4) (BB! J) ,J=1 *  NN ) 

115  MR  I TE  (  6 , 1 80C5 ) ( BB ( J) , J*1 , NN) 

122  MRI TE(6, 18001 1 

123  MS  =  MM  ♦  NN 

124  NS  =  MM  *  NN 

125  00  18110  I  =  1 , MM 

126  DO  18110  J  *  1,NN 

127  K  =  NN  *(I  -  1)  ♦  J 

1  30  CS(K)  =-CE  (  I  ,  J ) 

131  18110  AS! I , K )  =  0E( I  ,J) 

134  DO  18120  I  =  1 , NN 
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FORTH  AN  SOURCE  LIST  MATLIN 

ISN  SOURCE  STATEMENT 

125  DO  18120  J  =  l.MM 

1  36  I  I  =  MM  ♦  I 

127  JJ  =  NN  *  <J  -  1)  ♦  I 

1AU  18120  AS( I  I . JJ I  =  1. 

1  A3  00  18130  I  *  l, MM 

1A  A  18130  B S( 1 1  =  A A < I ) 

1A6  DO  1 81 AO  J  =  1 1 NN 

1A7  JJ  =  MM  ♦  J 

1*0  181A0  B  S(  J  J  )  =  B8 ( J) 

152  DO  18150  1  =  1  »  MS 

153  J  =  I  ♦  NS 

15A  IPATHSU)  *  J 

155  16150  AS(  I  ,  J  )  =  1. 

157  NNS  =  MS  ♦  NS 

C  -  SIMPLEX  ME THOO  SOLUTION  GF  LINEAR  PROGRAMMING  PKUBLEM 

1 6U  EPSLP  =  • 1 E-5 

C 

161  CALL  TODAY(0,I TIME , IUAT) 

C 

C  -  BEGIN  ITERATION 

C 

162  1820C  COS  T  =  0. 

163  DU  18205  1*1, MS 

It  A  18205  COST  =  COST  ♦  CTSU)*BSm 
C 

C  —  COMPUTE  ZS  AND  ZMCS  VECTORS 
C 

166  18216  00  1  822l>  J=  1  , NNS 

167  Z TS*0. 

170  DO  18217  1*1, MS 

171  18217  Z  T S*Z  T S  ♦  C TSI I ) *AS( I , J) 

173  ZS(JI=ZTS 

17A  18220  ZMCSI J»*ZS( JI-CSI JI 
C 

C  -  SELECT  MAXIMUM  ZMCS.  IF  NO  POSITIVE,  END. 

C 

176  18230  CMAX*ZMCSI1I 

177  JMA  X* 1 

200  DU  18250  J*2,NNS 

201  1F(CMAX-ZMCS( Jill  8240 ,182  50,182  50 

202  1 82 AO  CMAX=ZMCS(  J) 

203  JMA  X= J 
2CA  18250  CONTINUE 

206  IF( ZMCSI JMAXI-EPSLPI 18800,18800, 18260 
C 

C  -  SELECT  MINIMUM  DS ( I J =BS ( I ) /AS l I , JMAX )  WHERE  All, JMAX)  IS  PUSITIVt 

C 

207  ’,8260  DSMI  N*l.fc+35 

210  1 82 7G  DO  18350  1*1, MS 

211  IF(AS( I, JMAX) -.16-6)1 83 50, 18350, 18300 

212  18300  DST*BS(I)/AS(I ,JMAX» 

213  IF IOST-OSMIN 1  18310,183  50,18350 
21A  18310  OSMIN=UST 

215  I  MI N= I 

216  13350  CONI INUE 
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FORTRAN  SOURCE  LIST  MAT L IN 


,,ARMS»*L*MATLIN 

ISN  SOURCE  STATEMENT 


C 

C  -  CUMPUTt  Nth  MATRIX  ATS 

C 

*2u  DO  1840o  1=1, MS 

2  21  d  T SI  1  I = A  S ( I , JMAX )  *  BSIIMINI  /  AS (  IM I N, JMAX ) 

222  184CU  ATSI ( I ) = AS(  1  »JMAX)/AS( I  MIN, JMAX) 

224  TEMP- AS I I  MIN, JMAX ) 

225  TEMP2  =  ZMCS(JMAX)  /  TEMP 

226  THETA  =  BSIIMIN)  /  AS ( I  MI N, JMAX ) 

227  COST  =  COST  -  THETA  *  ZMCSIJMAX) 

23c  WRITE! 6,1 80 19ICQST 

231  00  18410  J  *  1 , NNS 

232  1 E41c  A  TSJ ( J I = AS  {  I  M I N , J I 

23*,  00  1 8526  1*1, MS 

235  IF( I-IMINI1  8450,18500,18450 

236  1  d4  50  IFIATSII II) 18455, 18525,18455 

237  1  E455  B SII) -B SI  1 l-BTS 1 1  > 

240  00  18475  J  «  1 , NNS 

241  IFIATSJIJII  18460, 184  75, 18460 

242  1846c  A  S<  I  ,  J )  =  AS(1,J)  -  ATSHI)  *  ATSJIJt 

2  43  IFIASI  I,J»-.lE-lul  18462,18462,18475 

244  18462  IFIASII.J)  -  .  IE- 10 )  18475 , 1846  5, 1  8465 

245  18465  A  S( I ,JI*0. 

246  1E475  CONTINUE 

250  GO  TO  18525 

251  ld50U  OU  18510  J-l.NNS 

252  AS! 1 , J)  =  AS(  I  , J) /TEMP 

2  53  IFIASI 1.JI-.  IE- 10) 18502,18502,1851C 

254  18502  IFIASI 1 , J )  -  .  IE- 10 >  1 8510 , 18505 , l B505 

255  ie505  ASI I , J)*0. 

256  1E510  CUNTINUE 

260  B S I  I  I *BS( I ) /TEMP 

261  18525  CUNTINUE 
C 

263  00  18550  J-l.NNS 

2t4  1E55U  ZMCSIJ)  «  ZMCSIJI  -  ATSJIJI  *TEMP2 

C  -  SUBSTITUTE  IPATH  OF  JMAX  FCR  1MIN,  C  UF  JMAX  Fuk  IM IN 

C 

266  IPATHSI  IN;?;; -JMAX 

267  f.TC{  I  MINI  «C  SI  JMm.v» 

C 

C - TRANSFER  BACK  TO  BEGIN  >  TE  RAT  I  ON 

C 

27C  GO  TO  1823G 

C 

271  1 88C0  CUNTINUE 

272  CALL  TUUAYll,!  TIME, I  OAT) 

273  TIME  «  FLOAT  1 1 T I  ME  ) /60  . 

274  WRl TEI6, 180151  TIME 

2  75  WRl TE I  6, 18020)1  IPATHSI I  I  ,  I «1 , MS > 

3C2  WRIT£l6,18005)IBSI II ,1-1, MS) 

30  7  WRITE  I  6,1 8605 1  I CTS I  I  1,1-1,  MS) 

314  CALL  EXIT 

315  ENO 
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it.  amtract  Thg  algorithm  described  in  this  paper  is  used  to  solve  a  special  class  of 
linear  programming  problems  characterized  by  constraint  coefficient  matrices 
having  generalized  transportation  structure.  Specifically,  n  available  resources 
are  allocated  to  m  capacity-limited  operations  (where  the  operational  capability  of 
assigning  the  jth  resource  to  the  ith  operation  is  known)  such  as  to  maximize  the 
total  profit  for  the  system.  The  row-and-column  structure  of  such  problems 
permits  an  algorithm  more  efficient  than  the  general  simplex  algorithm  to  be  used 
to  solve  moderate*  sized  problems  (problems  where  loop-tracing  techniques  or 
equivalent  schemes  are  not  required).  It  is  not  required  in  the  problem  statement 
that  all  the  resources  be  allocated  or  that  all  operations  be  performed  to  capacity 
limits.  It  is  characteristic  of  such  problems,  however,  that  the  optimizing  solu¬ 
tion  usually  requires  that  at  least  one  of  the  two  conditions  holds,  i.e.,  either 
supply  or  demand  is  exhausted.  The  paper  contains  a  description  of  the  algorithm, 
a  computer  program,  an  example  illustrating  its  application,  and  some  compari¬ 
sons  with  the  general  simplex  algorithm  in  solving  the  same  problem.  (  , 
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